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In Sec. 6.6, the text and equations following Eq. (6.107) until Eq. (6.111) should be 
corrected to read: 

...As pointed out in Ref. [44], Eq. (6.106) leads to an estimate for the Bose glass transition 
temperature from the high-temperature flux liquid phase, complementary to the estimate in 
Ref. [40]. We note that the average disorder potential for columnar defects is Vd = —U b 2 /d 2 , 
where U ss 0q/2(47tA) 2 denotes the defect potential well strength, b its width, and d the 
average distance between the columns [40]. Accordingly, the defect correlator becomes 
A = Uffi/d 2 . One important modification has to be taken into account near the local- 
ization transition (more precisely, for temperatures above the depinning temperature T^ p ), 
however, namely the effective thermal renormalization of the repulsive flux line interactions. 
In the boson picture, this can be obtained via a summation of the hitherto neglected vertex 
corrections, which in two dimensions leads to the replacement [11,51] 

V 

V ° ~" 1 + [W4tt(£;bT) 2 ] ln(l/nA 2 ) ' (6 ' 108) 
In the dilute limit n\ 2 <C 1 the scattering length a is therefore to be replaced with 

(k B T) 2 

a = —A > rr, x9 s • (6.109) 

47r ln(l/nA 2 ) 
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Upon introducing the characteristic temperature k^T* = b^iUo) 1 / 2 [40], using n ~ n 
B/(f) , and inserting Eq. (6.109) into Eq. (6.106), we find that 



7 W 



(t bg /tY 



(6.110) 



Thus C44 diverges as expected at the Bose glass transition temperature (compare Eqs. (D13) 
and (D6) in Ref. [40]), 

n l/4 

I I, . I ■ I 1 ; 

BG = J 



in qualitative agreement (apart from the logarithmic term) with a similar analysis in 
Ref. [44]. Although one should not trust the exponent for the vanishing of c^ -1 predicted 
by the approximate formula (6.110), this analysis suffices to estimate the transition temper- 
ature. ... 



The passage following Eq. (6.125), until Eq. (6.129), should be replaced with: 

... As for the linear defects, we can estimate a critical temperature from Eq. (6.124). The 
disorder correlator for planar defects is Ao = U^b 2 /d, where b denotes the diameter of the 
defect potential, and d the average defect distance. [Delete Eqs. (6.126) and (6.127).] Upon 
applying the replacement (6.109) again, Eq. (6.125) yields 



cl x ~\T) = {nh)- 1 1 - (Tb G /T) 



(6.128) 



and the instability occurs at the generalized Bose glass transition temperature marking the 
onset of localization to the defect planes 



T BG = T* 



In ^° 



3/2' 



m 2 d \ nB \ 2 b 



1/4 



(6.129) 



We are indebted to Johannes Kotzler for pointing out the errors in our previous estimates 
to us. 
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Abstract 

The influence of different types of disorder (both uncorrelated and correlated) 
on the superfluid properties of a weakly interacting or dilute Bose gas, as well 
as on the corresponding quantities for flux line liquids in high-temperature 
superconductors at low magnetic fields are reviewed, investigated and com- 
pared. We exploit the formal analogy between superfluid bosons and the 
statistical mechanics of directed lines, and explore the influence of the differ- 
ent "imaginary time" boundary conditions appropriate for a flux line liquid. 
For superfluids, we discuss the density and momentum correlations, the con- 
densate fraction, and the normal-fluid density as function of temperature for 
two- and three-dimensional systems subject to a space- and time-dependent 
random potential as well as conventional point-, line-, and plane-like defects. 
In the case of vortex liquids subject to point disorder, twin boundaries, screw 
dislocations, and various configurations of columnar damage tracks, we cal- 
culate the corresponding quantities, namely density and tilt correlations, the 
"boson" order parameter, and the tilt modulus. The finite-size corrections 
due to periodic vs. open "imaginary time" boundary conditions differ in in- 
teresting and important ways. Experimental implications for vortex lines are 
described briefly. 

PACS numbers: 05.30.Jp, 67.40.Yv, 74.60.Ge. 
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I. INTRODUCTION 

The theory of weakly interacting and/or dilute superfluid Bose systems has been studied 
for several decades, and now constitutes a well-established branch of condensed-matter 
physics (see, e.g., Refs. |]-|3||). More recently, upon utilizing the path integral representation 
of many-particle quantum mechanics |4]||, considerable progress has been reported for the 
even more difficult task to quantitatively understand the physical properties of strongly 
interacting boson superfluids such as Helium 4 below the Lambda phase transition line 
0. In addition to the detailed investigation of correlation functions in boson systems with 
these and other many-particle methods ]7|-|J, field-theoretic and renormalization-group 



approaches have been put forward in order to analyze these interesting systems |]T0Hl6f • 

While the statistical mechanics of superfluid bosons as such is thus a well-advanced field, 
historically there has been much less effort to study the influence of quenched disorder on 
superfluid properties; e.g., the effects of point disorder (even to lowest order in the defect 
concentration) on correlation functions. Characteristic quantities such as the condensate 
fraction and the superfluid density have been systematically evaluated only in the past few 



years |T7Hl9[ , in part motivated by the prospect of boson localization and related quantum 
phase transitions driven by the interplay of disorder and quantum fluctuations for superfluids 
in porous media pl|-P7f. 

Additional impetus for the investigation of disordered bosons stems from studies 
of magnetic flux lines in type-II high-temperature superconductors (for recent reviews, 



see Refs. pq , |29|| ), motivated by the formal analogy of the quantum mechanics of two- 
dimensional bosons with the statistical mechanics of (2+l)-dimensional directed lines 
|30|-p3[ . As is readily seen from the path integral formulation of the many-particle quantum 
mechanics for superfluids HH, the corresponding particle world lines behave as directed 
elastic strings with periodic boundary conditions in the direction of propagation z (imagi- 
nary time), with the boson mass m representing an effective line tension ei, and quantum 
fluctuations oc h mapping on thermal fluctuations oc T (see Sec. [IT]). This mathematical 
correspondence can be further exploited for the case of magnetic flux lines, as summarized 
in Table |, with the two-dimensional boson density n and conjugate chemical potential 
translating to the magnetic flux density B and external magnetic field H, respectively. 
Similarly, the (imaginary) velocity v and superfluid density p s correspond to a transverse 
magnetic field h and the inverse tilt modulus c±l of the flux line ensemble. 

As is depicted in Fig. [j], in a pure system, the superfluid phase of bosons is represented 
by a highly entangled line liquid R, while the normal-fluid phase corresponds to a disentan- 
gled liquid of flux lines. Finally, the hexagonal Abrikosov flux lattice is represented by a 



1 With the appropriate periodic boundary conditions, this reflects the occurrence of macroscopic 
particle exchange processes; compare the very suggestive pictures and discussions in Ref. II. 
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two-dimensional crystalline solid in the boson language []31|| . The prospects for finding an 
intermediate "supersolid" phase with both macroscopic phase coherence and crystalline or- 
der (not shown in the figure) may in fact be more promising for flux lines [p4| , |35| . If periodic 
boundary conditions are applied for the superconductor along the magnetic-field direction, 
the fact that in two dimensions a second-order Berezinskii-Kosterlitz-Thouless transition 
separates the normal and superfluid phases implies that there should be a sharp distinction 
between the entangled and disentangled line liquid states [the Kosterlitz-Thouless transi- 
tion temperature as function of the boson density T^(n) mapping on a critical sample 
length L(B) as function of the magnetic field]. However, if instead one assumes free rather 
than periodic boundary conditions, which should in fact be more realistic for most flux line 
problems, the correspondence between inverse temperature and sample size along the field 
direction is only approximate, and there need not be a sharp phase transition between an 
entangled and disentangled flux liquid | 36fl . This observation demonstrates that the issue 



of boundary conditions becomes important when the above correspondence is utilized, and 
that open boundary conditions may lead to very different physical behavior for L < oo from 
what would "naively" be inferred from the boson picture. 

For magnetic flux lines in high-temperature superconductors, disorder changes the phase 
diagram in essential ways. E.g., the presence of uncorrelated point defects destroys the 
crystalline order of the Abrikosov lattice [37], and may lead to a new disorder-dominated 



ground state, termed the vortex glass pq| . For the application of high-T c superconductors 
in a magnetic field, an effective mechanism for flux pinning is required, in order to minimize 
dissipative losses caused by the Lorentz-force induced vortex motion and thus retain a truly 
superconducting state, defined by a vanishing linear resistivity. For that purpose linear dam- 
age tracks have been produced by heavy-ion irradiation leading to columnar defects which 



turned out to be very effective in localizing the flux lines [25,25|. Theoretically, this corre- 
sponds to a two-dimensional boson localization transition; the new ground state of vortices 
pinned to columnar defects is called the Bose glass [|3"U|-[4"T[| (see also Ref. pEfl). We have 
recently further exploited the quantum analogy, in this case with disordered semiconductors 
(the Coulomb glass), in order to show that long-range intervortex repulsions in the Bose 
glass produce a highly correlated ground state and an additional reduction in vortex trans- 
port |K| . Such considerations may eventually allow controlled "tailoring" of defects in such 
a way that magnetic flux lines are bound most effectively. For example, the introduction 
of a certain amount of "splay", i.e., variation in columnar defect tilt angles was suggested 
to drastically improve flux pinning [H3-H3]. It was also demonstrated that columnar defects 



effectively increase the interlayer Josephson coupling in the pancake vortex regime [[|7j. In 
addition, the pronounced effect of planar disorder (as is frequently present in the high-T c 
ceramics in the form of twin boundaries) was studied ||48|| , and recently also related to the 
so-called peak effect |49j . The investigation of the effect of transverse magnetic fields trying 
to tilt the flux lines away from, say, an array of parallel columnar pins, leads to an interesting 



localization problem in non-Hermitian quantum mechanics [50 



In this article, we compare and contrast the behavior of boson superfluids and vortex 
liquids in the presence of various types of disorder, emphasizing effects due to the different 
boundary conditions. Coherent-state path integral techniques [f||| provide a unified frame- 
work to discuss both problems. The "open" boundary conditions appropriate to vortex 
liquids are illustrated in Fig. |TJ(c) , and differ from the periodic boundary conditions relevant 
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to real bosons by a set of independent integrations over the entry and exit points of the 



directed lines, see Sec. II B 



Some of the varieties of disorder treated here (to leading order in the disorder strength) 
are compared for vortex lines and superfluids in Table ||. Point disorder in three-dimensional 
vortex arrays corresponds to a rather unphysical random space- and time-dependent poten- 
tial applied to, say, superfluid Helium 4 films. However, columnar and planar pins in vortex 
liquids are directly relevant to boson superfluid films on disordered substrates. Columnar 
pins, of course, correspond to "point" disorder in two-dimensional superfluids, as would 
occur on an amorphous substrate, while random planar disorder for vortex liquids is rep- 
resented by random lines etched (say, via microlithography) onto an otherwise regular sub- 
strate. As part of this review, we rederive via the coherent-state method results for point 
disorder in superfluids |T7|JT8| , and produce new results for random line disorder perturbing 



two- and three-dimensional bosons, and random planes interacting with a three-dimensional 
superfluid. Among other results, we find that the superfluid density and condensate fraction 
in a zero-temperature superfluid film can be significantly reduced by depositing a set of lines 
with random positions and orientations in the substrate. 

To first order in the defect concentration, the influence of weak point disorder on corre- 
lation functions and thermodynamic properties of the flux liquid was calculated in the limit 
L — > oo by Nelson and Le Doussal [51,|52|]; in the present paper, we shall extend these previ- 



ous investigations, which are based on the path integral representation of the partition sum 
and the Gaussian approximation, to <i-dimensional generalizations of the disorder displayed 
in Table |J, as wen as two families of symmetrically tilted columnar defects and randomly 
splayed extended disorder. Our calculations of how disorder affects the superfluid density 
generalize a method developed by Hwa et al. to determine the effect of splayed columnar 
defects on the flux line tilt modulus [J44|j . Our investigation of the role played by differ- 
ent (namely periodic and free) boundary conditions (first discussed in Ref. |31[| ) elucidates 



the subtle differences in the underlying physics of boson superfluids and magnetic flux line 
liquids 0. We obtain explicit results for the finite-size corrections for a number of thermo- 
dynamic quantities, such as the boson condensate fraction n and the inverse tilt modulus 
c~ll (superfluid density p s ); these predictions should be subject to direct experimental tests. 
By invoking the boson analogy for a finite system "naively", i.e., using periodic boundary 
conditions, one would expect no and c~ll to decrease as a powers of 1/L, because thermal 
excitations — k^T ~ 1/L — reduce the Bose condensate fraction and the superfluid density 
at nonzero temperatures (see Sec. |111 Bp ; for open boundary conditions, however, it turns 



out that both these quantities increase (see Sec. |III C| )! This behavior arises because at the 
surfaces the probability distribution for the time-averaged position of a wiggling flux line 



is determined by the bosonic "wave function" |3T],[|2|, and is therefore larger than in the 
bulk, where this quantity is given by the "wave function" squared (a Gaussian for a free 
random walker); see Sec. [11 B| and Fig. |2|. Near the surfaces, the enhanced thermal wander- 



ing facilitates entanglement, and renders the system softer with respect to tilting, effectively 
increasing and the overall C44. 



2 In a different context, finite-size effects have also been studied by Marchetti and Nelson, using 
a hydrodynamic approach 
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Most of our final results on vortex liquids are summarized for weak, delta-function-like 
interactions between lines in the direction perpendicular to z. This restriction is easily 
relaxed by inserting the full Fourier-transformed pair potential [see Eq. (|2.22|) 1 into the ap- 
propriate formulas. However, our results for flux liquids are also limited to regimes where the 
interactions between lines are strictly local in z, the magnetic-field direction. Although this 
restriction facilitates comparisons with nonrelativistic bosons, it does limit straightforward 
applications of our results to high magnetic fields in superconductors. The most problem- 
atic quantity at high fields is the vortex tilt modulus C44, which has a large compressive 
contribution due to the interactions nonlocal in z. As discussed by Larkin and Vinokur ET 



for isotropic vortex liquids using a nonlocal generalization of the boson mapping the 
renormalized tilt modulus may be written as 

C U = — (l + -——) , (1.1) 



An \ 4ir\ 2 n s 

where A is the London penetration depth and n s denotes the superfluid number density. The 
first term is the compressive contribution due to nonlocality in z. Nonlocality is expected 
to be much less important in the second "vortex part" of the tilt modulus which is 
inversely proportional to the superfluid density of the equivalent boson system (]. We thus 



expect that our results may be used as an approximation to this second term of Eq. ( |1 . 1|) at 
high fields. The generalization of the Larkin- Vinokur formula to anisotropic superconductors 
with effective mass ratio m±/m z has recently been derived by Geshkenbein [j55|l , whose result 
we quote here for completeness 

This article is intended to be self-contained and pedagogical, and we therefore include 
sufficient introductory material. In the following section, we review the description of both 
the quantum mechanics of superfluid bosons and the statistical mechanics of magnetic flux 
lines by coherent-state path integrals. The important role of Galilean invariance (for bosons) 
and affine transformations (for flux lines) will be discussed. In Sec. fTJ we study pure sys- 
tems, first at zero temperature (i.e., L — > 00 for flux lines), and then also the lowest-order 
corrections for finite temperatures (finite sample thickness for the case of magnetic vortices). 



In Sec. [IV|, the disorder contributions to density, momentum, and vorticity correlations, con- 
densate fraction (depletion), and superfluid density are evaluated for bosonic systems (peri- 
odic boundary conditions) to lowest order in the defect concentration, and explicit formulas 
are presented for point, linear, and planar disorder in two and three space dimensions. In 
Sec. |V|, we discuss magnetic flux lines in high-temperature superconductors, for which open 
boundary conditions provide a better description, and determine the corresponding physical 



3 The correlation function which gives the superfluid density (see Sees. [II] and [III| below) depends 
only on fluctuations in the tangents to the vortex lines about configurations parallel to the z axis. 
Locality in z is an excellent approximation for these fluctuations, see App. A of the second part of 
Ref. g0|. 
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quantities, namely the magnitude of the "boson" order parameter, density and tilt correla- 
tions, and the renormalization of the tilt modulus by the disorder. Due to the generality of 
these computations the resulting formulas are rather lengthy; therefore we present explicit 
expressions for the cases of point disorder, splayed columnar defects, parallel extended and 
two families of symmetrically tilted defects in Sec. [VT]. Our findings for new kinds of cor- 
related disorder in boson superfluids, and the effects of different boundary conditions, are 
summarized and discussed in Sec. |V11| . In the Appendices, we derive the fluctuation formula 
of Hwa et al. [44] for the tilt modulus in very thick samples (in the thermodynamic limit) 
using affine transformations, discuss the "phase-only approximation" and its limitations, 
present a detailed derivation of the correlation functions and their disorder contributions for 
the case of open boundary conditions, and list some useful formulas required for performing 
the Matsubara frequency sums and momentum integrals. 



II. PATH INTEGRAL REPRESENTATION FOR BOSONS AND FLUX LINES 

WITH DISORDER 

In this section, we set the stage for our investigations of weakly interacting superfluid 
bosons and magnetic flux lines subject to disorder. By writing down the corresponding path 
integral representations for the partition sum, we elucidate the close formal relationship 
between the two distinct physical problems, and point out the differences. We also comment 
on the important role of Galilean invariance and affine transformations, respectively. 



A. Superfluid bosons 



Consider the Hamiltonian for iV bosonic particles in d space dimensions, with mass m, 
interacting via the scalar pair potential V(r), and subject to an "external" disorder potential 
Vb(r), 



n 



bos 
N 



N 

E 

i=l L 



— V 

2m 1 



Voir, 



1 



oEm^-r,!) 



(2.1) 



with the understanding that the corresponding iV-particle wave function is symmetric with 
respect to all possible permutations. As a consequence, the trace leading to the partition 
sum is to be restricted to these totally symmetric bosonic states {n} with energy eigenvalues 



E, 



bos 



^bos 
N 



E 



(2.2) 



bosonic states {n} 



where f3 = l/k^T. 

By applying the standard methods of rewriting the Hamiltonian (|2.1| ) in the language of 
second quantization and transforming to coherent states, the following imaginary-time path 
integral representation for the grand-canonical partition function is readily derived 0J|, 



7bos 



-s[Tp*,ip]/h 



(2.3) 
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with the action 

Sty*, ^ = f h dr I d d r Ur(r, r) ^lll + |1 | V ^(r, r)| 2 - L + Vb(r, r) 

+ - /rfVy(|r-r'|)|^(r',r)| 2 |^(r,r)| 2 



r,r 



(2.4) 



where /x is the chemical potential. The constraints for the fields ip and ip* in Eq. ( |2.3|) stem 
from the permutation symmetry requirement and thus reflect the boson statistics; they can 
be interpreted as periodic boundary conditions for an imaginary-time slab which at finite 
temperature is of "length" j3h. 

In the spirit of Landau-Ginzburg mean-field theory, we apply the method of steepest 
descent and find a nonlinear Schrodinger equation as the stationarity condition || 







5S[ip* 7 ip] _ dip(r 7 r 



n r i 

— V 2 ^(r, r) - n + Vb(r, r) ^(r, r) 
2m L J 



5^>*(r, r) dr 



(2.5) 



Upon specializing to short-range interactions with V^(r) = VoS(r), where Vq = V{<\ = 0) 
and V(q) = Jd d rV(r)e~ l(ir is the Fourier-transformed potential, taking the configurational 
average (denoted by an overbar), and assuming that |^(r, r)| 2 m n independent of r and r, 
this yields — (/i + Vd) + n Vo ^ l /no = 0, with solutions 



n = for // + Vb < 
«o = (/i + ^)/Vo for/i + VyJ>0. 



(2.6) 



The presence of disorder shifts the transition point to a phase with nonzero Bose condensate 
fraction no, which in a pure system would occur when u > 0. In the above mean-field 
approximation, the free energy in the superfluid phase becomes 



F gr (T, f2, u) = —k B T In Z, 



sir 



u + Vd Q/2Vq 



(2.7) 



where f2 = / d r denotes the volume of the system. Therefore the mean particle density is 



n = -n- L (dF gI /d/j,) 



n 



(2.8) 



in this approximation, and the mass density p ~ ran®. Similarly, the isothermal compress- 
ibility becomes 



Kt = n 2 (dn/dp) Tn » (njjV ) 1 , 
which implies for the macroscopic sound velocity (with ks ~ Kt) 

Cl = (ram^) -1 / 2 « (n Vo/m) 1/2 . 
Notice that ci vanishes both for no — * and Vq — > 0. 



(2.9) 



(2.10) 
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Eq. (|2.8|) implies that within mean-field theory, which neglects fluctuations, all particles 
enter the condensate in the superfluid phase. When density and phase fluctuations are 
taken into account, however, one finds that actually the condensate fraction is no < n (see 
Sec. |n| and Refs. []l],|2j,|4]]). The Bose condensate is formed by those particles which take 



part in macroscopic ring exchange processes on scales larger than the thermal de Broglie 
wavelength; the required symmetrization of the many-particle wavefunction then leads to 
phase coherence. In the path integral picture, superfluidity and tiq ^ can be related to 
a macroscopic fraction of world lines entangled and connected by virtue of the periodic 
boundary conditions, as is illustrated in Ref. 0. Note that n , the thermodynamic order 
parameter, is conceptually quite distinct from the superfluid density p s , which is a transport 
coefficient (see Sec. [Hi]) , defined as an appropriate limit of a dynamical response function 
in terms of the path integral representation, p s corresponds to the mean-square world 
line winding number ||. 

If we now allow for an inhomogeneous condensate [including an inhomogeneous chemical 
potential u(r, r)], and transform to real time according to r — > it, Eq. ( |2.5|) becomes 

It is then convenient to introduce the "polar" representation in terms of the new fields 7r(r, t) 
and 0(r, t), 



^(r,t) = ^n + n(r,t)e l@ ^ , (2.12) 

implying that the density fluctuations are described by the 7r field, 

n{r,t) = |^(r,t)| 2 = n + n(r,t) , (2.13) 

and furthermore assume small fluctuations from the mean condensate density only: ip fa 
y/n^(l +7r/2n )e ie [the average of the phase field Q(r,t) is taken to be zero]. We choose 
the constant n in Eq. ( [2.12|) to be exactly equal to the average boson order parameter 



squared, 

n = mr,t))\\ (2.14) 

which is a direct measure of the macroscopic occupation number of the p = momentum 
state. Note that in general n ^ (n(r,t)) = t)\ 2 ), although we will often set n ~ n, 

because (n(r,t)) w in the Gaussian approximation [see Eqs. ( |2.13| ) and (|2.^ )j. At T = 0, 
uq = n only for the ideal (non-interacting) Bose gas. 
Upon defining the superfluid velocity according to 

v,(r,t) = hVQ(r,t)/m , V x v s = , (2.15) 



the imaginary part of Eq. ( |2.11 ) then yields the continuity equation for the irrotational 
superfluid flow 0, 



dn(r,t) 
dt 



n(r,*)v,(r,*) =0 . (2.16) 
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On the other hand, from the real part of Eq. ([2.1 1|) the characteristic correlation length 

£ 2 = h 2 /4mn V = Ti 2 /Am 2 c\ (2.17) 

is inferred, and if the variation of the density 7r(r, t) is slow on a scale £, £ 2 V 2 7r <C it, the 
real part of Eq. (|2.11 ) becomes 



dv s (r,t) 



yv s (r,t) 2 



-V/i(r,t) (2.18) 



(note that it = 5n = r^K^b^ ~ Sfi/Vo), which is the fundamental equation of superfluid flow 
0. Finally, this may be cast into an alternative form, using the fact that superfluid flow is 
irrotational, V x v s = 0, 

<iv s (r,t) <9v,,(r,t) r . . . . . . . 

m ^ J = m s ^ } + [v.(r, t) ■ V] v s (r, t) = - V/i(r, t) , (2.19) 

which can be interpreted as Euler's equation for nonviscous flow, driven by a chemical 
potential gradient. 

In Sees. [TTT] and [VI we extend this analysis to incorporate both fluctuations in the fields 
7r and O (in the Gaussian approximation) and in the disorder potential Vo(r) (to lowest 
nontrivial order). First, however, we discuss the boson analogy for the statistical mechanics 
of magnetic flux lines in some detail. 



B. Magnetic flux lines 



In order to describe the large-scale properties of magnetic vortices in high-T c materials, 
which are strongly type-II superconductors, we work in the London limit, and define each 
flux line by its (two-dimensional) trajectory r^(z), as it traverses the superconducting sample 
of length L parallel to the external magnetic field H, which is taken to be aligned along 
the z direction. The ensuing Gibbs free energy for N such vortices subject to a disorder 
potential Vo(r), is [EH], 



G N (H) 




1 + 



m_i_ 



m, 



driiz) 



dz 



1/2 




rAz - 



^2.20) 



here, e\ = 6q h\K is the effective line tension, eo = (0o/4vrA) sets the energy scale, 0o = hc/2e 
is the magnetic flux quantum, k = \/£ is the ratio of the London penetration depth A and 
the correlation length £, and the material anisotropy is embodied in the effective mass ratio 
m±/m z (<^C 1 for high-temperature superconductors). Finally, the screened logarithmic 
intervortex repulsion is given by 



V{r) = 2e K (r/\) 



(2.21) 



where K (x) is a modified Bessel function, with the asymptotic behavior K (x) oc — hire for 
x — > 0, and K (x) oc x~ 1 ^ 2 e~ x for x — > oo. In Fourier space this potential reads 
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V(q) = V /(l + \ 2 q 2 ) , 



(2.22) 



where Vq = (fl/An. Note that we have assumed here that the vortex interaction is local in 
z, which is a good approximation when the flux lines are essentially straight and parallel 
to z, so that (\dr /dz\ 2 ) <C 1. In the same spirit, we can then expand the square root in 
Eq. ( 2.20Q , and reexpress the integral over the magnetic flux density b(r) in terms of 0o an d 
N, which leads to 



G N (H) = f iNL + F N [{r l (z)}] , 
where we have introduced the "chemical potential" [see Eq. fl2.26Q below] 



u = H ei= — (H 

P 4tt 4tt v 



H, 



(2.23) 



^2.24) 



which vanishes at the upper critical field H Cl = 47rei/0o, and the vortex free-energy func- 
tional 



Fs[{Ti(z)}] 



N 



i=l 



drAz) 



dz 



D 



TdZ) 



-J2v(Hz)-r,(z)\] 



(2.25) 



with the modified line tension ei = {m±/m z )ei. 

A full statistical treatment requires a summation of e~ GN ^ kBT over all possible vortex 
trajectories {r^z)}; thus the partition sum becomes 



7 R _ ST 1 „!iNL/k B T 7 ft 



N=0 



N ) 



(2.26) 



where the canonical partition function for a system of N lines reads 

N 



Z N = f[ /P[r l (z)]e-^^«>]/ fcBT . 
i=i 



[2.27) 



With Eq. ( |2.25| ) this already has the form of a quantum-mechanical partition function in 
the path integral representation, i.e., for the world lines of N particles of "mass" e\ moving 
through "imaginary time" z, interacting with potential V(r), and subject to a disorder 
potential Vd(t). According to Eq. ( p.24| ), the external magnetic field plays the role of a 
chemical potential, and at T = vortices will start to penetrate the sample when H > H Cl ; 
furthermore, the areal particle density is related to the magnetic flux via n = B/fo. 

To render this quantum analogy of the classical statistical mechanics problem of directed 
lines more precise []3T| , |52] , we note that the partition sum ( 2.27 ) can be rewritten in terms of 
the transfer matrix e~ H N L / k ^ T connecting neighboring constant-z slices, i.e., as an integral 
over a quantum-mechanical matrix element (in "imaginary time"), 



N . N . 

Z% = II \K II / ■ ■ ■ r' N \e-^W\ ri ...r N ) 

i=i J i=i J 



[2.28) 



where the states |rj • • ■ r N ) and (r^ ■ • • r^| describe the entry and exit points of the vortices, 
respectively, and the Hamiltonian is [compare Eq. (|2.1|) 1 
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N 

^jv = E 

i=i 



i^V(|r ( - rj |). (2.29) 



By inserting a complete set of many-particle energy eigenstates |n) in Eq. ( |2.28| ), and intro- 
ducing the zero-momentum state 

N 

| Pl = 0--- PA r = 0) = n Jdr i \T 1 ---T N ) , (2.30) 

i=l 



we find 



^= (Pi = • • • p N = 0|e-^ L / feBT | Pl = • • • P7V = 0) 



^|(n|pi = 0---p N = 0)\ 2 e -< L ' h * T . (2.31) 



The "quantum-mechanical" interpretation of this matrix element is that of a system pre- 
pared in the ground state |pi = ■ ■ ■ pjv = 0) of the ideal Bose gas at z — 0. It is then 
allowed to evolve in imaginary time under the influence of interactions and disorder. After 
"time" L, this interacting superfluid is projected back onto the ideal Bose gas ground state. 

To see that only boson states are involved in the statistics of these fictitious quantum 
particles with "free" boundary conditions, note that the zero-momentum state fl2.30| ) is 



symmetric under permutations, and therefore only those eigenstates of the permutation- 
symmetric Hamiltonian ( |2.29| ) contribute to the partition sum ( |2.31| ) that are permutation- 
symmetric themselves. Hence the summation effectively includes bosonic many-particle 
states only, 

Z%= E |^|pi = 0---p^ = 0)| 2 e-^° s W. (2.32) 

bosonic states {n} 



Upon comparing with the "pure" bosonic partition function ( |2.2j ), we see that there is 
a difference arising from the weights involved in the projection onto the zero-momentum 
ground state. Yet in the limit L — > oo, only the lowest-energy bosonic state contributes to 
the sum ( |2.32|) , and upon identification of the respective quantities in Table |, the quantum 
mechanics of superfluid bosons at T = becomes completely equivalent to the statistical 
mechanics of directed lines in the thermodynamic limit. If we impose periodic boundary 
conditions on the vortices, i.e. use a toroidal geometry in the field direction, the exact 
correspondence of Eqs. ( |2.32| ) and ( |2.2| ) even survives for finite sample lengths (flux lines) 
/ nonzero temperatures (bosons), and a coherent-state path integral representation for the 
grand-canonical partition function ( |2.26| ) may be derived, which is precisely of the form 
of Eqs. ( p.3|) , (|2.4j) . Note that thermal fluctuations oc k B T correspond to quantum fluctua- 
tions oc h, while the torus length L maps onto the inverse boson temperature (5h. In this 
case, a sharp Berezinskii-Kosterlitz-Thouless transition is to be expected from an entangled 
(superfluid) phase at large L to a disentangled (normal fluid) phase at small L pO|j3T| . 

For magnetic flux lines in most superconducting samples, however, periodic boundary 
conditions are rather artificial, and a more adequate description of finite systems may be 
obtained using open boundary conditions as embodied in Eq. ( |2.32| ), where the line exit and 



entry points are integrated over freely, corresponding to an ideal Bose gas with 
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V>(r, 0) = ip(r, L) = y/n 



(2.33) 



Note that this rigid constraint at the boundaries means that the "depletion" n — uq = 
(\ip\ 2 ) — \ (ip)\ 2 (see Sec. |IIIB| ) vanishes identically at the top and bottom of the sample. As 
is shown in Ref. jnj], within the mean-field approximation the sole effect of the projection 
amplitude in (|2.32|) is a shift of the chemical potential according to 



k B T . N 



(2.34) 



where Z\ represents the partition function for a single flux line. Physically, this amounts 
to a downward entropic renormalization of the upper critical field H Cl [see Eq. ( |2.24| )], due 
to thermal fluctuation of vortices. With the replacement /i — > jl, and the correspondences 
of Table |, the results of the previous Sec. |II A| can then be taken over. Fluctuation effects 
may furthermore be described by a coherent-state path integral with an action of the form 
(|2.4|), provided the additional constraints on the fields if) and ■0*, or n and 0, respectively, 
are explicitly taken into account, see Ref. plj and Sec. [Ill C| below. In this geometry, the 



exact correspondence between L and (3h in a real boson system is lost, and the sharp phase 
transition from an entangled to a disentangled flux liquid may in fact be replaced by a 
smooth crossover |3C| . We also remark that the bosonic field operators ^ and i\> correspond 



to magnetic monopole and anti-monopole creation operators in the flux line picture. 

The vortex free energy ( |2.25| ) was derived for a constant magnetic field H in the z 
direction. The response to "tilt", i.e., to a small slowly varying magnetic field Hj_(r, z) 
perpendicular to z is also of interest. Such perturbations couple directly to the "tangent 
field" of the flux lines, 



tr,z 



y, dv.i z) 

i=l 



dz 



5 m 



(2.35) 



and the free energy (|2.25| ) is modified according to [] 



Iq 



~[ Jo dz 



(2.36) 



where 



v(r,^) = — -H ± (r, z) 
47rei 



(2.37) 



4 To derive this result, note that the perpendicular component of the magnetic field h±(r,z) 
in Eq. fl2.20p satisfies the London equation, b_|_(r,2i) = (m 2 /mj_)A 2 V5_bj_(r, z) + A 2 <9fb_|_(r, z) + 
4>0 t(r, z). Upon substituting for b^(r, z) in Eq. ( 2.20 ), the derivative terms on the right-hand side 
can be eliminated via an integration by parts with negligible error, provided the variations in the 
external field Hj_(r, z) are sufficiently slow. For a treatment of the response to more general field 
perturbations, see Appendix A of Ref. |5l|. 
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The field v(r, z) thus couples to the flux lines like a vector potential in this "Lagrangian" 
formulation of the physics. In the transfer matrix representation of the vortex partition 
function, the Hamiltonian ( |2.29| ) becomes |51 



N 



N 



knT. 



N 



1 



+ £Vb(r i ) + =5>(|r,-r i |) 

i=i z i^j 



(2.38) 



Note that the tilt field generates an imaginary vector potential in this Hamiltonian formula- 
tion. If Hi is constant, then tilting the external field acts like an imaginary Galilean boost 
with velocity iv = i(fi H±/4:7iii on the fictitious bosons. 

As a simple, but useful application of the mapping of vortex physics onto quantum 
mechanics, we briefly consider the thermal wandering of a single flux line (N = 1) in a 
sample of length L, starting at the origin r = 0, and ending at the opposite surface at 
r = rj_ |]40|,|4"2l. With Eqs. (|2.28|) and ( |2.29| ), the constrained partition function for this 



situation reads 

Z(r ± , 0; L) = (r ± \e- n i L/k * T \0) , (2.39) 
and the probability distribution for the vortex tip position at the upper surface becomes 



V(r ± ; L) = Z(r ± , 0; L) / j d d r ± Z(r ± , 0; L) 





0) (r ± 


n) 


e -E^L/k B T 




0) Jd d r ± (r ± 


^ e -E%L/k B T 



[2A0) 



where in the second equation we have inserted a complete set of energy eigenstates |n) with 
eigenvalues E^. For thick samples, L —>■ oo, the ground state dominates the above sums, 
and therefore 



V(rr,L) 



V>o(r±) 



Jd d r ± ip (r±) 



-(Ef-E*)L/k B T\ 



(2.41) 



here, V'o( r ±) = ( r ±|0) is th e ground state "wave function" with eigenvalue Eq, and Ef is 
the energy of the first excited state. The probability distribution for the flux line to have 
wandered a distance r± at the surface is thus proportional to the ground state wave function 
[note that ^o( r ±) is positive and nodeless]. 

On the other hand, if we are interested in the more general problem of one flux line 
entering the sample at r = rj and leaving it at Tf, and ask what is it its average thermal 
wandering probability distribution V(r; L) at height z far away from the two boundaries, 
the answer is 



V(r;L) = Z(r;L) / Jd d rZ(r;L) 



(2.42) 



where, using the definition ( |2.39| ), 

Z(r; L) = J d d n J d d r f Z(r f , r;L — z) Z(r, r i5 z) . (2.43) 
Upon inserting a complete set of energy eigenstates, for thick samples the final result is 
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V(r; L) = Jf° (r)2 ^ [l + O U'^-^W)] , (2.44) 

i.e., in the bulk, this probability distribution is proportional to the square of the ground state 
wave function. Note that for periodic boundary conditions along the z direction, which rule 



out any special surface effects, Eq. ( 2.44|) describes the thermal wandering in the entire 
sample. On the other hand, for open boundary conditions the surface behavior is different 
from that in the bulk, and Eq. Q2.41|) applies. We therefore conclude that near the surfaces, 
thermal flux line wandering is enhanced for open boundary conditions, as depicted in Fig. |2|. 

We comment finally on the meaning of the "boson" order parameter for directed lines. 
On a formal level, superfluidity is usually accompanied by long-range order in the correlation 
function 



specifically 



G(r,z ] r',z') = ^(r,z)r(r',z')) , (2.45) 



lim G(r,z;r',z') = (^(r, z)) <^*(r', z 1 )) = n + , (2.46) 

|r— r'|^oo 



where the average is evaluated via coherent-state path integrals. To understand what this 



long-range order means for vortices [|52[], it is helpful to consider first a situation where 
the order is manifestly short-ranged, in the Abrikosov flux lattice. Figure ||| shows that the 
operators i\> and ^ create vacancy and interstitial strings in the vortex crystal, terminating in 
"magnetic monopoles" . The composite operator in Eq. Q2.45| ) creates an extra line at (r', z') 
(i.e., a column of interstitials in the solid), and destroys an existing line at (r, z), creating a 
column of vacancies. The lowest-energy configuration is then a line of vacancies (for z' > z) 
or interstitials (for z' < z) connecting the two points with an energy as proportional to the 



length s of this "string". It follows that the correlation function ( 2.45Q decays exponentially 



to zero (i.e., oc e -°1 r - r 'IABT^ f or i ar g e separations in the crystalline phase. In an entangled 
line liquid, on the other hand, vacancies and interstitials "melt" away and the asymptotic 
string tension a vanishes, implying long-range order (|2.46|) in G(r, z; r', z'). In this limit, 
{ip(r,z)) = (ip*(v',z')) = e ~ Einon / kBT , where E mon is the energy of an isolated "magnetic 
monopole" . 



C. Galilean invariance of the pure system and affine transformations 

Before we proceed with the investigation of fluctuation effects in boson superfluids and 
flux line liquids, we discuss the fundamental symmetries of the pure systems, namely in- 
variance with respect to Galilean transformations and uniform tilts, respectively. We begin 
with the disorder-free boson Hamiltonian Q2.1|) , and study the effect of a Galilean boost 
with constant velocity v = const., 

r' = r + vt , t — > t' = t , 
P' = P + Nm\ , 

E' = E + v • P + l -Nmv 2 , (2.47) 
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where P is the total momentum, and E the total energy of the iV-particle system. 

In the coherent-state path integral representation (with imaginary time r = it), the par- 
ticle density n(r, t), momentum density g(r,r), particle current j(r, r), and energy density 
H(r, t) = H(r, r) — jin{r, r) read 



n(r,r)=|^(r,r)| 2 , (2.48) 

^*(r,r)V^(r,r) -^(r,T)Vf (r,r) , (2.49) 

h 2 ,_., ,, , i 



v r,rj= mj(r,r) = — 



^ (r ' T)= 2^ IW(r ' r)|2 "^ ( ^ (2 ' 50) 



A Galilean transformation ( |2.47| ) is then represented by a unitary transformation acting on 
the fields, 

</>(r, r) -> V'(r', r') = ^(r, r) e ^ v - r A ; (2.51) 

as is easily seen, this leaves the particle density invariant, 

n'{r', t') = n(r, r) , (2.52) 

while the transformed momentum and energy densities become 

g / (r',r')= g(r, r) + mvn(r,r) , (2.53) 
W'(r', r')= W(r, r) + v • g(r, r) + irrn; 2 n(r, r) , (2.54) 



in accord with Eqs. (|2.47|) . Note that the third term in Eq. (|2.54 ) may be absorbed into a 
shift of the chemical potential, 

H'(r / ,r / )=H(r,r)+vg(r,r) , (2.55) 

//= li - mv 2 /2 . (2.56) 

In the polar representation 

^(r, r) = Jn + n(r,r) e te ^ , (2.57) 



the effect of the Galilean transformation is 



n'(r',r')=n(v,r) , (2.58) 
e'(r',r')= 0(r,r) + mv • r/H . (2.59) 



Upon introducing the superfluid velocity according to Eq. ( [2.15|) 



v s (r,r) = liV9(r,r)/m, g(r,r) = mv s (r,r)fi(r,r) , (2.60) 

one finds 

v;(r',r')=v s (r,r)+v. (2.61) 
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For the case of directed line liquids, the Galilean boost ( |2.47 ) translates into an affine 
transformation, namely a uniform tilt away from the z axis, 

, , , dr dr' dr h . 

r -> r ' = r + hz e x , z -»• z' = z , - > — = — + —, (2.62) 

az az az e\ 

with h = const., which may be interpreted as a transverse external magnetic field. As dis- 
cussed above, the uniform tilt h/e^ = 0oHj_/47rei corresponds to a Galilean transformation 
with imaginary velocity iv (see Table |I|). This imaginary boost velocity means that the 
phase does not change in the simple fashion described by Eq. fl2.59| ). In the corresponding 
quantum-mechanical problem one has to consider the effect of imaginary vector potentials 



leading to the appearance of non-Hermitian operators ||50|| . The transformation formulas for 
the energy density may be read off from the previous expressions; the shift in the chemical 
potential becomes 

// = // + h 2 /2i 1 . (2.63) 



We finally remark that under a Galilean boost ( j2-47|) or affine transformation ( |2.62 ), the 
disorder potential is modified according to 

Vb(r(r), r) -> V D (r' + ivr', r') , (2.64) 

and 

V D (r(z), z) - V D (r' - hz'/^, z') , (2.65) 

respectively. Therefore, disorder breaks "Galilean" invariance; however, for uncorrelated de- 
fects for which the disorder correlator is a constant, the symmetry is recovered by taking the 
quenched disorder average. In the Appendix |A], we review how these affine transformations 
can be exploited for the calculation of the disorder contribution to the line tilt modulus in 



the thermodynamic limit L — > oo ||44j| . A method which is more convenient for finite L is 
discussed in the next section. 



III. FLUCTUATIONS AND CORRELATION FUNCTIONS IN PURE SYSTEMS 



We are now ready to evaluate correlation functions for our model action (|2.4| ). We 
shall henceforth regularly use the boson notation and language, but nevertheless keep the 
application to flux liquids in mind as well. We start by rederiving and summarizing the 
known results for the pure system with periodic (bosonic) boundary conditions 0, before 
we extend our analysis to open boundary conditions (Sec. [Ill (J| ), which are more appropriate 
to flux line systems, in the presence of uncorrelated and correlated disorder (Sees. |TV|,[V|). 

A. Correlation functions in the Gaussian ensemble with periodic boundary 
conditions (Bogoliubov approximation) 

To separate density and phase fluctuations, it is convenient to use the polar representation 
(p.57 ) for the fields ip and ip*. Note that the path integral in the partition sum ( |2.3|) is 
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originally defined by integration over Reijj = + it cos and Imip = ^n + it sin B; 
the Jacobian corresponding to the above transformation is 1/2, and being constant may be 
absorbed into the path integral measure. As we shall only be dealing with small fluctuations, 
we extend the integration range for the phase field from [— 7r; 7r] to [— oo;+oo]. Thus 
the grand-canonical partition function becomes a functional integral over the independent 
variables it and 0, 



7bos 



n{r,/3h) = ir(r, 0) 
e(r,[3h) = 0(r, 0) 



D[7r(r,r)]X?[0(r,T)]e 



-(S [7r,e]+S int [7r,e])/ft 



(3.1) 



where, upon splitting the disorder potential into its mean value and the fluctuations around 
it, Vb(r, r) = Vd + SVd(t,t) and using Eq. ( |2.6| ) in the superfluid phase, /i + Vd = n Vo, 
the harmonic and "interacting" parts of the action read 



S [w,O] 



dr / d rl ikir^, r) 



de(r,r) h 2 



Or 



[V7r(r,r)f 



h 2 no 

2m 



[V0(r,r) 



SV D (r, r) 
~n 2 Q V (3h{l 







n + 7r(r,r) 


AS 





d d r' V(\r-r' 



it r , r 7r r, T 



and 
Swt[n, 0] 



dr d d r< 



imno n + 7r(r, r) 



[Vvr(r,T)f 



— 7r(r,r) [V0(r,r)f 



(3.2) 



(3.3) 



respectively. This separation of the action amounts essentially to a splitting into single- 
and multi-quasiparticle excitations. In this paper, we shall entirely neglect the anharmonic 
terms (|3.3|) , and restrict ourselves to the Gaussian approximation. Note that taking into 
account the nonlinear fluctuation terms ( |3.3| ) perturbatively would amount to an expansion 
in powers of h 2 (~ T 2 for flux lines); this would then include further contributions from the 
disorder potential which enter via their linear coupling to the density fluctuations 7r. 

The harmonic action ( p.2|) is of course immediately diagonalized via (discrete) Fourier 
transformation, 



1 



with the inverse 



f H f 0K 

V(q,w m ) = JdrJ Q d,Ti/){T,T)e 



-i(q-r— oJ m r) 



(3.4) 



(3.5) 



In the thermodynamic limit, the discrete wavevector sum is as usual replaced by an integral, 
J2q . • • Q(2ir)~ d J d d q . . .; the periodic boundary conditions in imaginary time on the other 
hand enforce the discrete bosonic Matsubara frequencies to be 

2tt 



m 



(3h 



m 



0,±1,±2,... 



(3.6) 
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which, similarly to the wavenumbers, approach a continuum as (3 — > oo, J2m ■ ■ ■ ~~ > 
f3h/(2n)- 1 JZ du.... 

For pure systems, i.e., if 6Vc(r, r) = 0, the effective action in the Gaussian approximation 
can be written in the matrix form 



n z v phn 



5o ^ e] = ^ q 5J 0( " q '" Wm) ^ ( " q '"^ m) ) A(q ^ m) (°(q,Cj) ~^ ' 

(3.7) 

with 

M ,_(n h 2 q 2 /m -huj m \ , , 

A(q, u m ) - y ^ v{q) + h 2 q2/4mno ) ■ (3-8) 

The density and phase correlation functions can therefore be obtained by simple inversion 
of the matrix (|3.8|) , 

* A -if n , . \ _ ( me B (q) 2 /n h 3 q 2 uj m \ . . 

here, e-Q^q) denotes the Bogoliubov quasi-particle spectrum 



e B (q) = ^n h 2 q 2 V(q)/m + (h 2 q 2 /2m) 2 , (3.10) 

which at long wavelengths displays a linear dispersion, 

6 b(<?) - > hciq for q — > , (3-H) 

where c\ is the speed of sound from Eq. ( |2.10| ). 

We can assign a physical meaning to the Bogoliubov spectrum by investigating the 
density-density correlation function, 

S{r, r; r', t') = (n(r, r)n(r', r')> - n 2 « (7r(r, t)tx{y\ t')) . (3.12) 

Note that by adding a source term to the action, i.e., by performing a Legendre transforma- 
tion according to S[8fi] = S[tp*,ifj] + ^ h dr Jd d r 5/x(r, r)|^(r, r)| 2 , we can also write this in 
the form of a functional derivative, 



s(r,T-yy) = h 2 6 lnZ - m 



a %(r,r)> 5 „ 



6[Sfi(r,T)]S[S t i{T',T')] 



(3.13) 

5)1=0 



which shows that the imaginary-time correlation function is actually a quantum-mechanical 
response function. From Eq. ( |3.9|) , in Fourier space we immediately find as a consequence 
of translation invariance in space and time 

) = S{q,u m ) (27i) d 5(q + q!)PhS m ,. m/ , (3.14) 

with the reduced density correlation function 
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n hq 2 /m 
u& + e B (g) 2 /^ 



g(q,^ m ) = - - 2 ; (3.15) 



by performing the simple Matsubara frequency sum [Eq. ( p2| ) in App. p 1| | we can readily 
transform back to imaginary time, 

1 r>nh 2 n 2 p(P- T / h ) € v{q) _i_ p e B (?)T7i 

with the static limit 

5(q '° ) = 2^) COth ^' (3 - 17) 

which is precisely the static structure factor in the Bogoliubov approximation for weakly 
interacting Bose gases |],|2|]. Note that S(q, 0) is inversely proportional to the excitation 
spectrum when (3 — > oo. If we now analytically continue to t = —it, and carry out the usual 
Fourier integral, we arrive at 

S(q,u,) = / dtS(q,t)e*« = [e*BM*( W - e B ( g )/h) + *( w + e B (g)/ft 



(3.18) 



i.e., the dynamic structure factor directly measured in scattering experiments. One there- 
fore expects quasiparticle peaks with the dispersion relation ( |3.10| ), thermally excited with 
a Bose-Einstein distribution (note the detailed-balance factor e@ e ^ q '). The Gaussian ap- 
proximation for the density and phase fluctuations produces exactly the same results as 
the familiar Bogoliubov approach in many-particle quantum mechanics, where a canonical 
transformation for the creation and annihilation operators is employed in order to diagonal- 
ize the second-quantized quasi-particle Hamiltonian 0,0]. In the long-wavelength limit, the 
elementary excitations of the Bose condensate are therefore no longer of the free-particle 
form €o(q) = h q 2 /2m, but display the linear dispersion ( |3.11| ) and may thus be interpreted 
as phonon modes. Depending on the form of the interaction potential V(q), the disper- 
sion relation may display the roton minimum characteristic of strongly interacting Bose 
liquids like superfluid Helium 4. However, our harmonic approximation is only valid either 
for weakly interacting systems where V(q) is sufficiently small, or for dilute Bose gases for 
which the pair potential may be strong, but where the scattering length a = mVo/An is 
small compared to the average interparticle spacing. In that case, one can exactly sum all 
the ladder diagrams describing multiple scattering events, with the result that eventually 
mV(q) is to be replaced by the effective s-wave scattering amplitude f = mV the 
dimensionless coupling in the Bogoliubov theory, serving as an expansion parameter in the 
weak-coupling limit, is in fact g 2 oc n^ 2 ^ d Vo fll |. In both situations the quasiparticle 



dispersion is a monotonically increasing function, and shows no roton minimum. 

We also remark that a second route to arrive at the result fl3.18| ), is to remember that 
S(q,u m ) = hxn(q,^ m ) is the density response function; by analytic continuation in fre- 
quency space, oj m — > i{uj +irj), r\ — > + , and splitting x n (q, oS) = X^(q, to) + ix'ni^ a; )' where 
Xn(q^) = [Xn(q,^ + ^) + Xn(q,^-^)]/2 and x'n(<l,u) = [xn(<liUJ + ir))-Xn(<liU-ir))]/2i 
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denote its components which are continuous and discontinuous across the real axis, respec- 
tively, one may then infer the dynamic structure factor (in real time), i.e., the corresponding 
correlation function, from the quantum-mechanical fluctuation-dissipation theorem 

2hZ(q,uj) = (l - e~^) S(q, u>) . (3.19) 

As a second important correlation function, we consider the momentum (mass current) 
correlator in the Gaussian approximation, defined as 

C*,(r,r;r',r') = (^(r, r)g 3 (r', r')> - (ft(r, r))(^(r', r')> . (3.20) 



According to Eq. ( |2.60| ), g = mj = Ti{riQ + 7r)VB; inserting this into Eq. ( |3.2U| ), we can 
factorize the four-point correlators into all possible two-point functions (i.e., apply Wick's 
theorem), while the three-point correlators vanish in the Gaussian ensemble. In Fourier 
space, this yields 

C -(q,w TO ;q',£«O = -ft 2 nQg^(e(q, w TO )e(q', u' m )) 

PiPj (7r(q- P,^m-n)7r(q' - p / ,^ m '_ n '))(6(p,^ n )6(p / ,u;„/)) 



+ (?r(q- p,u; m _ n )e(p',^„/))(e(p,u;„)7r(q' - p',a; m /_ n /)) 



(3.21) 



Note that this quantity is not simply proportional to the phase correlation function as 
obtained from Eq. (|3.9|), 



(e(q,c m )e(q',^)) = m T^ 2 { n fS (27r)^(q + q0^^ . (3.22) 

[In the Appendix [B], we discuss the "phase-only approximation" in some detail, which 
neglects the additional correlators in Eq. (|3.21|) 1. Using the Gaussian correlation functions 
to be read off from Eq. (|3.9| ), one finds for the transverse current response function 



X±(q,w ro ) = — J 1 — g Pg'(q)C^(q, u m ) 



1 ^ gV-(qp) 2 (q~P) 2 e B (p) 2 A +p 2 Um-nUn ,„ „„, 

~ (d - l)/?n P V„ q 2 P 2 [<-n + e B (|q - pl) 2 /^ 2 + e B ( P ) 2 /h 2 ] ' { ' ' 

Here, the transverse projector is -P^(q) = Sy — P^(q), with P«(q) = <Mj/<? 2 - Note that 
Xj_(q, Mm) would have vanished entirely, had we not taken the nonlinear terms in Eq. ([3.21 



into account, because the contribution from the phase fluctuations in Eq. fl3.21| ) is purely 
longitudinal Q. 



5 The above factorization of the four-point correlation function actually goes beyond the Gaussian 
approximation (tree level). In field-theory language, it corresponds to a computation of the one- 
loop contributions to the composite-operator cumulant ([7rVj6][7rVjG]). Therefore, in order to 
consistently evaluate the longitudinal current response function to this order, one would have to 
include the one-loop corrections to the phase correlations as well, which would involve the nonlinear 
vertices of Eq. (3.3). In the limits u> m — > and q — > 0, the loop terms actually cancel, and the 
ensuing result for the total density is p = lim^o Xy (q, 0) as n^m. 
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In the limit u m — > (first), followed by q — > 0, the transverse current response func- 
tion yields the norma-fluid density, for this quantity is defined precisely as the transport 
coefficient relating to the system's response to transverse motion This yields 



Pn = limx±(q,0) 



(3.24) 



upon performing the Matsubara frequency sum [Eq. (|D12| ) in App. p 1| , we find 

,2 „Ad' ^ 



Pn 



(3h A f d d q 



Ad J (27r) d Vsinh[/5e B (g)/2] 



(3.25) 



and therefore p n (T = 0) = 0; at finite "temperatures", Eq. (|3.25|) will be evaluated explicitly 
in two and three dimensions in Sec. 1111 B. 



It is interesting to note that the transverse momentum response function is actually inti- 
mately related to the vorticity correlation function, which (in d = 2 and d = 3 dimensions) 
is defined as 



(3.26) 



V i:j (r, r; r, t') = m 2 ( V x g(r, r) . V x g(r', r) 



In Fourier space, and after factorizing according to Wick's theorem, this becomes in terms 
of the density and phase fluctuations 



Vij(q, uj m ; q',uj m >) — eiki£jk>i>- ( — ^ ^ 
k,i,k',i' \rnp\t) Ptil j n p / ;UJn , 



Yl (?* _ Pk)Pi(q'k> - P'M> 



x 



(7r(q - p, u m - n )7i(q - p', u m > (0(p, u n )Q(p', uj n >)) 
+(7r(q- p,a; m _ n )e(p / ,u; n /))(e(p,u; n )7r(q / - p',u m ,_ n ,)) 



(3.27) 



Following the same procedure as above, one may demonstrate explicitly that the vorticity 
correlations are purely transverse, Vy(q, u m ) = Pij(q)Vij(q,u m ) = 0, and after some 
algebra finds 



m 2 V ± ( y q,u m ) = hq 2 x±(q,uJm) ■ 



(3.28) 



It is therefore the intrinsic vorticity fluctuations that render the transverse current response 
function and the normal-fluid density Q3.25| ) nonzero at positive temperatures. These fluctu- 
ations would be confined to the vortex cores in a "phase-only" Kosterlitz-Thouless picture 
of finite-temperature superfluid films. 

Finally, it is natural for flux lines to introduce the "tilt" correlation function by consid- 
ering an action perturbed by an imaginary "velocity" field v(r, r) given by Eq. ( |2.37| ), 

r/3h r 

S[v] = S[ip*,ip] + J drjd d r [iv(r,r) -g(r,r) - mv(r,r) 2 n(r,r)/2] (3.29) 
[compare Eq. ( |2.54j) 1, and then defining [5lJ 
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T- 



fr,r; r', r') 



h 2 5 2 lnZ gr [v] 



m 2 5vi(r,r)Svj(r',r') 



(t i (r,r)* j (r',T')) 



v=0 



(nh/m) Sij8(r — r')5(r — r') — Cy(r, r; r', r')/m 2 



which in Fourier space becomes 

7ij(q,a; m ) 



(nh/m) Sij - Cij(q,u m )/m 2 . 



(3.30) 



(3.31) 



Thus, tangent-tangent correlations for flux lines are intimately related to momentum corre- 
lations in the equivalent superfluid. To extract the tilt modulus from Eq. ( p. 31 ), one has to 
take q — > first, followed by u m — *> 0. These limits mimic the actual physical situation for 
open boundary conditions, reflecting the response to an external tilt field which acts on the 
top and bottom faces (perpendicular to z) of a superconducting slab p3| , [51~]l . We can again 
isolate a matrix of "transport" coefficients, namely 



n h) lim Ty(0, uj r , 



1 

nm 



5i L_ 

13 nm(3VL 



E PiPj 



6 B ( P f/h 2 



e B (pf/h 2 + ujl 



(3.32) 



For an isotropic system this matrix reduces to the tilt modulus, c"- 1 = c v u 1 5 i j, with 

(nm) 2 (nm - p n ) _1 = p 2 /p s , (3.33) 



-li 



where p s = p — p n is the superfluid density, and we have set p = nm. In the flux liquid 
notation, p n = for (disorder-free) thick samples (L — > oo), and thus p s = p = nm and 
C44 = nm = nei in the vortex language. 

At this point we mention that in an Abrikosov flux lattive the transverse tilt correlations 



assume the form 51 



Tj_ (q, uj„ 



n 2 hu) 2 m 



CmT + nmuo^ 



(3.34) 



which reduces to the transverse part of the flux liquid result [cf. Eq. ( |3.31| )1 when the shear 
modulus c 66 vanishes. Thus, in the crystalline phase, c v ^ = nm as in the liquid, but p s = as 
the result of the appearance of the shear modulus and the opposite order of limits (uj m — > 
to be taken first, then q — > 0). 

As we shall discuss in detail in Sees. [TV] and |V|, the effect of disorder will be twofold. 
First, for disorder correlated in the xy plane, the tilt modulus tensor (|3.32|) will in general 
be anisotropic (i.e., the normal and superfluid densities, being transport coefficients, will 
depend on the direction of motion with respect to the correlated defects); second, disorder 
will lead to a renormalization of p n and hence p s and c\y A phase transition into a localized 
phase of bosons / flux lines will occur when eventually the superfluid density and therefore 
the inverse tilt modulus vanish (in a certain direction at least). This "generalized Bose 
glass" phase then constitutes a true superconductor with zero linear resistivity P0 , while 



the flux liquid with p s 7^ behaves like a normal metal |j28 |. 
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B. Green's functions, condensate depletion, and normal fluid density for a superfluid 

Bose gas 

It is useful to define the following Green's functions for the Bose system (with periodic 
boundary conditions in imaginary time), 



G(r,T-r',r')= (i>{r,T)r{r',T')) -n 
G 12 (r,r;r',r')=^(r,r)^(r',r'))-n 



(3.35) 
(3.36) 



(compare Refs. PH). For a translationally invariant system one deduces from Dyson's 
equations the general structure 



G(q,u m ) = iu m + hq 2 /2m - fi/h + X n (-q, -u m ) /D(q,u ri 



(3.37) 



which defines the self-energy £n(q, w m ), while the anomalous Green's function Gi 2 (q, oj Ti 
takes the form 



G 12 {q,u, 
and the denominator reads 

1 



-Si 2 (q, uj m )/D(q, uj t , 



(3.38) 



D(q, LO r , 



+ 



lUJ r , 



~ ^( S n(q>^m) - Sn(-q, 



-Ur. 



hq 2 /2m - jj/h+ -(£n(q, u m ) + £n(-q, -ay, 



Si 2 (q,o;) 



(3.39) 



For noninteracting bosons, £u(q, u m ) = S i2 (q, 0J m ) = 0, and hence the anomalous Green's 
function vanishes, while the free propagator becomes 



G (q, to m ) 1 = -iu m + hq 2 /2m - fi/h . 



(3.40) 



Furthermore, it can be shown that S n (0, 0) — S 12 (0, 0) = fi/h to all orders in perturbation 
theory |L|. 

In order to calculate these Green's functions in the harmonic approximation, we expand 
Eq. ( |2.57| ) to first order in the density and phase fluctuations, 



^(r,r) 



-i , nr,r) , . 7r(r,r) 

1 + — + 26(r,r) - 

2n 



2n, 



-7r(r,r)0(r,r) - -Q(r,rf 



(3.41) 



To leading order, this yields for the Green's functions in Fourier space 
G(q 



1 i 
- — (7r(q,u; m )7r(q / ,u m >)) - -<7r(q, UJ m )e(q', u m >)) 
4n 2 



Gi 2 (q, uj m ; q', ay, 



+-(6(q, u m )ix(q\ u m i)) + no(0(q, o> m )9(q ow)) , 

1 i 
- — (7r(q,o> m )7r(q',aw)) + -(7r(q, u m )Q(q', u m >)) 
An 2 

i 

+~(Q(q,o> TO )7r(q / ,aw)) - n (O(q, a m )6(q', ow)) . 



(3.42) 



(3.43) 
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Upon inserting the results ( |3.9|) from the previous section, one readily finds for the (reduced) 
Green's functions in Fourier space [see Eq. Q3.14Q 1 



iu m + hq 2 /2m + n V(q)/h 
co 2 m + e B (q) 2 /h 2 



G(q,w m )= — : j-^ —2 , (3.44) 



G 12 {q,u m )= -— 2 , (3.45) 

u 2 m + e B {q) 2 /h 

and therefore with /x = noVo 

Su(q,w m )= «o[Vb + V{q)]/% , Ei 2 (q,u; m ) = n V(q)/h , (3.46) 
D(q,^ m )=^ + e B (g) 2 A 2 • (3-47) 

Here e B (q) denotes the Bogoliubov spectrum, already introduced in Eq. (|3.10D . The results 
( p.44| )- (|3.47| ), obtained in the Gaussian approximation, are again identical to those from 
the Bogoliubov theory applicable to either weakly interacting or dilute Bose systems IHpf. 
Another useful representation for the Green's functions consists in their separation into 
single poles and their residues, 

G(q.^)= . "5 M/f + - f' 2 M/f . (3-48) 
G 12 (q,u m )=-u(q)v(q)(— 1 p— + , 1 — ) , (3.49) 

with the weight functions 

«(g) 2 1 = 1 / n y(g) + ft 2 g 2 /2m 
^) 2 J 2 I e B (g) 



±1 , (3.50) 



obeying u(g) 2 — v (q) 2 = 1 and u(q)v(q) = noV(q)/2e B (q)] these are precisely the coefficients 
of the Bogoliubov transformation from the original boson to quasiparticle creation and 
annihilation operators. 

We can now use the Green's function to calculate the average particle number in the 
grand-canonical ensemble; carefully taking the normal ordering of the original Bose field 
operators into account ||||, we have to evaluate (with rj I 0) 



N - N = J d d r ((^(r, r - r/)^*(r, t + T]))- n ) 

= /tS* WW ^ e-^^VS^^VG^ Um . _ q? Wm/) . (3.51) 



nun 



(2vr) d (^) 

If in addition time translation invariance holds, we find 

(note the crucial exponential factor stemming from the time ordering). Assuming that 
the particle density n is known, Eq. (|3.52j ) can be regarded as a formula for the squared 
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magnitude of the boson order parameter, n = | {-^(r, t))| 2 . Upon using Eq. ( |3.48|) and 
Eq. (|D1| ) from App. |D 1| , we find the condensate depletion in terms of the quasiparticle 
weights ( |3.50| ), 



n — n 



d d q 
(2tt)° 



v(q) 2 + 



u(q) 2 + v(qf 

P /3e B (<?) _ 1 



(3.53) 



Because of the relation /z = n V , this actually constitutes an implicit equation for the 
chemical potential. The first term in Eq. ( 3.53 ) represents the zero-temperature depletion, 
i.e., the fraction of particles that are pushed out of the condensate as a consequence of inter- 
actions and quantum fluctuations; the second contribution yields the additional depletion 
caused by thermal excitations Q. 

For short-range interactions, V(q) ~ Vo = Ana/m, where a denotes the (s-wave) 
scattering length, and Eq. ( ^.53[ ) can be evaluated explicitly; at T = 0, we have (with 
x = hq/ ^8Tm a) 



n(T = 0)-n = J 



d d q 



v(qf 



(2ir) d 

POO 

(2an /h 2 ) d/2 T{d/2r 1 / x d - 2 



(2 + x 2 ) 1 / 2 -x-(2 + x 2 )- 1/2 



dx . 



(3.54) 



In d = 1, this integral diverges logarithmically at the lower limit, which means there is no 
Bose condensation in one dimension, in accordance with the Hohenberg-Mermin-Wagner 
theorem. For d = 2, the integral is equal to 1/2, and hence 



d = 2 : n(T = 0) — n = n a/h 2 , 
while in three dimensions |l|||,(|] 

d = 3 : n{T = 0) - n = (8/3Vn)(n a/h 2 f /2 



(3.55) 



(3.56) 



Note that Eq. (|3.55| ) with Ti — > k^T, m — » e± and Vo = (fro/^ 71 yields the "boson order 
parameter" squared for vortex lines as function of temperature in the limit of infinitely thick 
samples (L — > oo), 



no(T) =n 1 + e! {<p /Aixk B TY 



-i 



(3.57) 



The reduction in this quantity from unity is a measure of the confining effects and diminished 
entanglement induced by the repulsive vortex interactions. 



6 The depletion can also be derived somewhat more directly by using the definition ( |2.14| ): for, 
Eq. (|3~4l| ) implies that n = |(V>)| 2 ~ + (tt) - (vr 2 )/4n + i((ir9) - (6tt))/2 - n (@ 2 ), and 
hence because of n = (\tp\ 2 ) = uq + (it) one finds n — no = (vr) = (7r 2 )/4no — i((7T0) — (07r))/2 + 
no(@ 2 ) [compare Eq. ( |3,42j )l, which immediately leads to Eq. ( |3.53| ). Note again that formally the 
contributions to these composite-operator averages amount to loop integrals. 
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Returning to the boson representation, we now calculate the lowest -order finite- 
temperature corrections explicitly, approximating the Bogoliubov spectrum by its phonon 
branch, eB(qO ~ ficiq. Then 1 + 2v(q) 2 m mci/hq, and n(T) — n = An(T), with 



An(T) - I d q 1 + 2v ^ ~ r(d-l)C(d-l)m d _, 

where Eq. ( pi3| ) from App. |D 2j was used. In two dimensions, the ( function diverges, 
and hence a nonzero boson order parameter is only possible at zero temperature; at finite 
temperatures, there is still a superfluid phase characterized by a finite superfluid density p s 
(see below), and the phase transition will be of the Berezinskii-Kosterlitz-Thouless type. 
For d = 3, Eq. ( HTSgP yields HgJi 

d = 3 : An(T) = m(k B T) 2 j\2Tc'c 1 . (3.59) 

We can also utilize the Green's function (|3.35|) to rederive the normal-fluid density 
( p.25| ) following a famous argument due to Landau [0,0,0. Consider a superfluid confined in 
a cylindrical pipe, whose walls are moved along its axis at constant velocity v with respect to 
the superfluid. In this situation, the superfluid velocity is zero, and the normal fluid velocity 
is equal to v, because according to the two-fluid model only the normal fluid component is 
dragged along by the moving walls. Thus 

(g(r, t)) = m(j(r, r)) = p n v ; (3.60) 

the task now is to calculate the average current to first order in v, which then, by comparison 
with Eq. ( p.6C| ), yields p n . The definition ( |2.49| ) leads, for spatially translation-invariant 
systems, to (with rj j 0) 

(g(r ' T)) = / (S^W ^ e-^ + ^ e -^G(q, u m] -q, u m ,)_ v , (3.61) 



rn.rn' 



which reduces further to 

if translation invariance in imaginary time holds as well, which is of course the case for a 
pure system with periodic boundary conditions. The index "— v" serves as a reminder that 
we need the Green's function for a system with nonzero relative velocity between the walls 
and the superfluid. 

We can readily obtain the Green's function for moving walls from the known quan- 
tity G(q, u) m )o for fixed walls, Eqs. ( |3.44|) , (|3.48|) , by performing a Galilean transformation 
with constant velocity —v. The Matsubara frequencies in the transformed system are [see 
Eqs. 



tfV, r')= V(r + «vr, r) = -L £ ^(q, u m )e^-^-^ = £ ^'(q', k/J^^ 

^ q^m l'i w 'm 

q -»• q' = q > w m ^^ = o; m -iq-v; (3.63) 
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notice that the periodic boundary conditions are also obeyed by the transformed fields, 
ip'{r' -i-v/3h,/3h) =ip'(v',0). 

Upon inserting the shifted Matsubara frequencies in Eq. ( |3.48| ), we find 

rt \ l + v(q) 2 v(q) 2 . * 

G(q, w TO )_ v = — ■ rT7 r + ■ . w , 3.64 

-iu m + e B (q)/h - q • v iu rn + e B (q)/h + q • v 

with the unchanged Bogoliubov coefficient v(q), Eq. ( |3.50|) . Performing the Matsubara 
frequency sum [Eq. (pl|) 1 finally yields Landau's formula 

[the terms oc v(q) 2 vanish, because their combination is antisymmetric under q — > — q]. 
Eq. (|3.65|) can be interpreted as a sum over the momenta times a Bose-Einstein quasiparticle 
occupation number in a system with moving walls (where the energy eigenvalues are reduced 
by ftq • v). An expansion to linear order in v then again leads to the result ( |3.25| ). 

We now evaluate this expression in the phonon approximation, e B (q) ~ hciq. Via inte- 
gration by parts, the integral in Eq. ( |3.25| ) becomes of the form ( pi3j ), and hence 

pn[1) 2 d ~H d / 2 dT{d/2) h d ci +2 ' 1 ' 

Thus, in two dimensions we find 

d = 2 : p n = 3C(3)(k B Tf/2ixn 2 c\ , (3.67) 

with C(3) w 1.202, while for d = 3 HgQ 

d = 3 : p n = 2tx 2 (k B Tf/Ahh z c\ . (3.68) 

For a flux liquid with periodic boundary conditions in the z direction along the magnetic 
field, the leading finite-size corrections to the tilt modulus thus are (d = 2) 

Ac^(L) = -3((3)k B T/27r(ni 1 ) 2 ciL 3 . (3.69) 

Thus, C44 increases and it is harder to tilt vortices in finite-size systems with periodic bound- 
ary conditions, as one might expect. In the following section, we determine the correspond- 
ing finite-size corrections for the condensate fraction no and the tilt modulus (normal-fluid 
density) for the more realistic case of open boundary conditions. 

We finally remark that a crucial point in the calculation of the normal-fluid density is 
that both phase and density (amplitude) fluctuations are taken into account properly. It is 
tempting, given that the phase fluctuations are the dominant low-energy modes, to integrate 
out the 7r fields in the harmonic action ( |3.2|) , and then derive thermodynamic quantities like 
the depletion and the normal-fluid density from the ensuing effective action. As is already 
apparent from Eq. (|3.23| ), however, essential parts of the transverse current correlations will 
then be missing, e.g., there will be no vorticity in the system, unless vortices are explicitly 
introduced by hand. Hence one finds p n = even for T > in this approximation. This 
unsatisfactory result arises because it is actually not the leading terms (in a long-wavelength 
expansion), but the next-to-leading contributions containing the amplitude fluctuations, 
which produce a non-zero vorticity and normal fluid density at finite temperatures. The 
phase-only approximation (for the disordered boson system) and its limitations are more 
extensively discussed in App. |B[ 
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C. Finite size corrections for open boundary conditions 



We now proceed to study finite-size corrections for the flux line problem with free bound- 
ary conditions. To facilitate comparison with the periodic boundary conditions in the previ- 
ous section, we retain the boson notation. According to Eq. Q2.33Q , there can be no density 
and phase fluctuations at the bounding surfaces for all r, 

7r(r, 0) = 7r(r, ph) = = 6(r, 0) = 6(r, /3h) , (3.70) 

reflecting "ideal Bose gas" boundary conditions at the top and bottom of the sample. 
Eq. |3.70| amounts to a constraint for the path integrals as in Eq. (|2.3|) that has to be 



taken into account explicitly. Upon introducing the vector notation 



T(r,r)= ( tjlZl) , (3-71) 



7r r, r 



the two-point correlation functions with the constraints ( p.70|) requires evaluating the path 
integral 

JV[T(r, r)] T(r, r)T T (r', r') e~ s ™ ft. <T(r, 0) 
(T(r,r)T (r ,r ))= 



JV[T(v, r)] e-sm/n n r <f(T(r,0) y 
JV[\(r)]V[?(r, t)] T(r, r)T T (r', r') e -S[T]/R+i/«f r A(r)T(r) 



JP[A(r)]P[T(r, r)] e -5m/a+i/^A(r)T( 



r) 



(3.72) 



where in the second line a static auxiliary field A(r) has been introduced to enforce the 
constraints. In the Gaussian approximation, the correlation functions obeying the boundary 
conditions (|3.70|) can then be reduced in a straightforward manner to those with periodic 
boundary conditions that are encoded in the matrix A _1 (q, u m ) of Eq. ( |3.9|) . For the pure 



case, this has been worked out in detail in Ref. |]3T|, with the result 

(T(q, u m )T T (q , u m >)) = (2n) d 5(q + q')/iA _1 (q, u m ) f3h6 m) -. m > - A(q)A~ 1 (q, -u m , 

(3.73) 

where A(q) denotes the inverse of A _1 (q, r = 0) = (jSh)^ 1 J2 m A _1 (q, w ra ), 

fn hY/me B ( q ) 2 ^2^tanh^M . (3.74) 

VM; y me B {q) / n h q 2 J 2 v ; 

In App. 0, the analogous formula is derived for a general disordered system, and the pure 
result follows from Eq. ( |C9| ) by simply setting the disorder correlator to zero. Note that 
translation invariance only holds in the d perpendicular directions, while the second term 
in Eq. ( |3.73| ), representing the corrections from the open-boundary-conditions constraint, 
explicitly breaks translational invariance along the direction of the lines. 

Upon multiplying the matrices ( |3.74j ) and ( |3.9j ), we find for the density correlation func- 
tion with open boundary conditions 
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S(q,« m ;q / l w m O = (2n) d 6(q + q') — ^^-^ x 



x 



u4 + e B (g) 2 A 2 
2 ^ tanh ^ eB '- 9 '' WmU;m ' + e B(<?)V^ 2 



e B(g) 



2 ^, + 6 B ( ? ) 2 A 2 



(3.75) 



The frequency sums required for the equal-time correlation function follow from Eqs. (P5| 
and (|D8| ) , with the result 



S(q,r;q',r) 



E 



-i{Um+U m l)T 



S(q, uj m ; q', u m .) = S(q, 0) (2vr) d <5(q + q') , (3.76) 



%,0) 



n h 2 q 2 /3e B (g) 



2me B (q) 



cosh 2 [/3e B (g)/2]/ 2me B (q) 



n h 2 q 2 taBh @ e M 



(3.77) 



Compared to the corresponding result for periodic boundary conditions fl3.17| ), the hyper- 
bolic cotangent function has been replaced by a hyperbolic tangent, as if the excitations on 
the superffuid ground state were governed by a fermionic rather than a bosonic distribution 
function [cf. Eq. flD2f) ]; accordingly the limit 5"(q) = n h 2 'q 2 j '2me B (q) for f3% — > oo is ap- 
proached from below for open boundary conditions, as opposed to from above for periodic 
boundary conditions. 

Similarly, the Green's function (|3.42|) reads 



G(q, u m ; q', w m ,) = (27r) d 5(q + q) x 
hq 2 /2m + n V(q)/h 



eB(q) 2 /h 2 



^ + c B (?) 2 /^ 



iuj m (3h5 ri 



2h l3e B {q) uJmMW + e B (g) 2 /fr 2 ' 

^ 2 <, + e B (g) 2 /n 2 , 



es(g) 

2e B (g) 



ft 



tanh 



2 ^ + e B (g) 2 M 2 



(3.78) 



from which we can infer the average particle number with Eq. ( 3.51|) using the Matsubara 
frequency sums in App. p 1| ; instead of Eq. ( |3.53| ) for periodic boundary conditions, the 
result for the depletion now is 



n — n - 



d d q 
{2n) d 

d d q 
(27r) d 



v(q) 2 + [u(q) 2 + v(q) 2 
u(q) 2 + v(q) 2 



1 



ePe B ( q ) _ i s inh[ ( 5e B (g)] 



v ( q y - 



e Pe B (q) _|_ l 



(3.79) 



i.e., the finite- "temperature" (finite-size) corrections change sign, and the corresponding 
"excitations" are governed by a Fermi-Dirac distribution (with zero "chemical potential") 
rather than a Bose-Einstein function as in the case of periodic boundary conditions. 

With the aid of Eq. (P14j) , we can evaluate the "depletion", which is now actually a 
condensate enhancement, in the phonon approximation. Remarkably, in two dimensions the 
open boundary conditions lead to a finite correction, 



d = 2 : An(T) = -(In 2)m(k B T) / 2nh 2 



(3.80) 



in contrast to the logarithmic divergence one obtains for the case of periodic boundary 
conditions. The possibility of long-range order at finite temperatures in two dimensions 
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arises, due to the static surface "field" A(r) which enforces the open-boundary-conditions 
constraints at r = and r = (3%. This field strengthens the order parameter at the free 
surfaces, which promotes condensation in the bulk, when the thickness in the imaginary-time 
direction is finite. For higher dimensions we find 

d>3. An(T)-^l--j 2,-1^/2^/2) ^-2 i ( 3 - 81 ) 
thus, in three dimensions, 

d = 3 : An(T) = -m(k B T) 2 /24h 3 Cl , (3.82) 

which is (—1/2) times the depletion for periodic boundary conditions, Eq. ( |3.59|) . As we 
have already discussed in the Introduction, the physical reason for the signs in Eqs. fl3.80| )- 
( j3.82j ) is in the flux line language that thermal wandering is stronger at the surfaces, and 
thus line entanglement is facilitated, see Fig. ^. 

For the same reason, with open boundary conditions one expects the superfluid density 
p s (and the inverse tilt modulus C44) to be larger at finite (3fi than in the "zero-temperature" 
limit, and thus the finite- "temperature" / finite-size corrections should have opposite signs 
compared to the bosonic case with periodic boundary conditions [Eqs. (|3.25| ) and (|3.66|) 1- By 
repeating Landau's argument for tilted lines, we obtain the Green's function from Eq. ( |3.78| ) 
upon replacing uj m — > uj m — iq ■ v and u) m i — > u m i + iq ■ v [cf. Eq. (|3.63|) 1 . Then expanding to 
first order in q • v and performing the Matsubara frequency sums with the aid of App. |D 1| , 
yields the normal-fluid density, which has now become a function of imaginary time, 

ft ( smh[((3h-2T)e B ( q )/fi) _ 0h \ 

Pn[T) ~ d J (2ttY q \ sinh[/5e B (g)] 4 sinh 2 [/56 B (g)/2] ) " 

In the center of the bulk (r = (3h/2) this reduces to 

Mt ^ V2 )._§/^(_^) 2 , ,3.84) 

which is just Eq. ( |3.25| ) with a negative sign. This negative result does not of course have 
the usual interpretation of two-fluid superfluid hydrodynamics, i.e., of a physical "particle 
density" , once flux line boundary conditions are applied. In the phonon approximation, we 
find explicitly 

r(d + 2)c(d+i) (k B T) d+1 

Pn[ } ~ 2 d -^dT(d/2) h d ci +2 ' 1 ' 

and in two and three dimensions one has 

d = 2 : p n = -3C(3)(A;BT) 3 /27rft 2 c!f , (3.86) 

d = 3 : p n = -2tt 2 (k B T)y45h 3 cl . (3.87) 

For the finite-size correction of the tilt modulus with open boundary conditions this 
implies (in the flux line notation, d = 2) 

Ac^(L) = 3C(3)£;BT/27r(nei) 2 c?L 3 . (3.88) 

As anticipated, the system with open boundary conditions can be more easily tilted than 
one with periodic boundary conditions [Eq. ( p.69|) 1. 
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IV. DISORDER CONTRIBUTIONS FOR SUPERFLUID BOSONS (PERIODIC 

BOUNDARY BONDITIONS) 



In the preceding section, we studied defect-free systems of superfluid bosons and flux 
liquids, and thereby established notations, introduced the relevant quantities, explained our 
approximations, and discussed the relevance of boundary conditions. We now determine 
the influence of weak disorder on density and phase correlations, the depletion, and the 
normal-fluid density (tilt modulus). We start by investigating weakly interacting or dilute 
superfluid Bose gases, or flux liquids with periodic boundary conditions along the magnetic- 
field direction. The case of open boundary conditions more appropriate to vortices will be 
treated in Sec. [V|. We shall derive our results for general disorder first, specializing at the 
end to specific impurity correlators (Sec. [VQ) . 



A. Density, current, and vorticity correlations 



We retain our approximation to work in the Gaussian ensemble described by the action 
only, neglecting the nonlinear terms Q3.3| ). In Fourier space, the harmonic action, 
including the disorder terms, reads 



5oK9]= 



1 



E 



(6(-q, -U) m ) , vr(-q, -U m )) A(q, u m ) ( W " 



2(6(-q, -u m ) , 7r(-q, -u n 





6V D (q, u} r , 



- -n 2 V (3fin , (4.1) 



where the harmonic coupling matrix A(q, w m ) is given by Eq. fl3.8|) . In order to obtain the 
two-point correlation functions to first order in the disorder correlator 



A(q,u; m ;q / ,u; m /) = 5V D (q, UJ m )5V D (q', u m >) , 
we introduce the shifted fields (note that the Jacobian of this transformation is 1) 



(4.2) 



/ Q(q,u m ) 
\ 7r(q,a; m ) 



6(q,cu m ) 
7r(q,u; m ) 

0(q,u; m ) 



A 1 (q,w m ) 
6V D (q,uj n 





5V D (q,u m ) 
( u m /% 



" ^7r(q,a; ro ) j u) 2 m + e B (q) 2 /h 2 \ n q 2 /m 
upon performing the quenched disorder average we find 7 



(4.3) 



(6(q,u; m )0(q / ,u; m O) : 



me B {q) 2 /n Q h 3 q 2 
u4 + e B (g) 2 A 2 



(2Tr) d 5(q + q')Ph5 mj 

w m uj m i /h 2 



< + eB(q) 2 /hV m/ +e B (q'y/h 2 ^ 



A(q, u m ; q',uj r , 



(4.4) 



7 For a derivation of these results using the replica trick, see Ref. [51], App. B 
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(0(q, w m )ir(q',u m >))-- 



(27r) d <5(q + q') / 3M m! _ 



ul + e B (qy/h 2 ^ 
n Q q' 2 uj m /hm 

R + ^(g) 2 A 2 ][^ + e B (g') 2 A 2 ] 



A(q,^ m ; q,u r , 



(4.5) 



(7r(q, u m )7r(q',u m >))-- 



n hq 2 /m 



n^q 2 q' 2 /m 2 



(27r) d <5(q+q')/3M m ,_ 



A(q,u m ; q',u r< 



(4.6) 



Eq. ( |4.6| ), of course, gives the dynamic structure factor in harmonic approximation to 
first order in the defect correlator. In most situations, translational invariance in space and 
time will be restored statistically, i.e., 



A(q, oj m ; q', u m i) = A(q, u m ) (27r) d <5(q + q') (3M r , 



(4.7) 



in that case, one finds that the disorder contributions to the structure factor have the typical 
Lorentzian-squared form, 



S{q,u ri 



riQ%q 2 /m 



+ 



noq 2 /m 



ul + e B (qy/h 2 ' \ul + e B (qy/h 2 



A(q,u m ) 



(4.8) 



The corresponding static structure factor becomes, using Eqs. ( p5[ ) and ( |D9| ) of App. p 1| . 



— — n h 2 q 2 (3e B {q) 



n q- 



E 



A(q,u; n ) 



m J ph^[ul + e B {qY/h 2 Y 



(4.9) 



With the help of Eqs. ( |3.21| ),( |3~2"7D and some tedious algebra, we can also compute the 
current and vorticity correlations to first order in the defect correlator. Upon assuming that 
the system is statistically translation-invariant and isotropic, we can define longitudinal 
and transverse components of the current correlation function, and as in the pure case the 
transverse current response function is closely related with the [purely transverse) vorticity 
correlations, 



m V±(q,u m )= kq x_j_(q,a; fl 
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(q - p) 2 ^ 2 



P UJ m — n UJ n 



uJ 2 +e B (qy/h 2 



A(p,u n ) 



(4.10) 



For general, "time-dependent" configurations of disorder, the average mass current does not 
necessarily vanish; e.g., uniformly moving defects "drag" a certain fraction of the particles 
along, even at T = 0. Indeed, upon using Eq. (|4.5|), one finds 
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(q, uj m ))= — J2 P( 7r (q - P> 

' p,UJn 



.np 1 



E 



n )0(p,cu n )) 
p(q-p) 2 w n A(q-p,w„ 



6 B (|q- P |) 2 /^]K + eB(p) 2 A 2 



(4.11) 



Upon returning to Eq. ( |4.1(J| ), we find for the normal fluid density in an isotropic system 



limv±(q,0)= - , 
^° ^p,^ [e B (p) 2 /n 2 + 



Ef 2 



6 B (p) 2 A 2 -^ 2 



n 



dmf3TiVt 



V- 4 eB(p) 2 A 2 -3w 2 
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e B (p) 2 A 2 + ^ 2 



(4.12) 



we shall evaluate this for interesting special cases in Sec. [IV (J . For static disorder, i.e., 

A(q,w m ) = A(q)/?M m , , (4.13) 
Eq. ( |4.12| ) further simplifies to 



_ 

Pn 



d d q 



Ad J (2vr) d Vsinh[^e B (g)/2] 



n 
dm 



d d q ( hq 



(2n) d le B (g), 



A(q) • 



(4.14) 



An identical formula for the disorder-induced renormalization of the normal-fluid density 
was derived by Giorgini, Pitaevskii, and Stringari | 18fl . Note that the defect contribution to 
Q4.14T ) is temperature-independent, see Sec. [VB . 



For flux liquids, Eqs. ( |3.31| ) and ( |4.10| ) lead to the vortex contribution to the tilt modulus 



„v -i 



n 2 h) 1 lim T i:j (0,u m )= [S. 



nm 
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6 B (p)7^ 2 + ^ 2 



E PiPjP 2 
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A(p,u r , 



(4.15) 



which for an isotropic system agrees with Eq. (|3.33| ). We also remark that the last term in 
Eq. ( [4.15|) can formally be written as 



QiQjQ 
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d d q 



E 



e B (q) 2 /h 2 -3uj 2 



d d q 

nr J (2vr) d 

2n~^E I (27p m Ph ¥ A(q ' dul 



[e B (qy/h 2 + ul 
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A(q,w„ 



d 2 



S(q, u m ) 



(4.16) 



where S(q,cj m ) is the density correlation function of the pure system ( |3.15|) . Eq. ( [4.15| ) is 
thus the generalization to finite systems of the tilt modulus formula of Hwa et al. |[44|| . In 
App. |A], we review how the defect contribution to the tilt modulus can be obtained using an 
affine transformation; it is, however, not straightforward to use the same idea and method 
to obtain finite-size corrections. In the following subsection, a more general expression for 
the normal fluid density and the tilt modulus will be derived, valid even when translational 
invariance along the r direction is not restored statistically. 
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B. Condensate fraction and normal fluid density 



Eqs. ( j4.4j )-( |4Tq ) yield the disorder-averaged Green's function in harmonic approximation 
( P-42j ), which reads, to first order in the disorder correlator, 

— r iuJ m + hq 2 /2m + n V(q)/h d 

G(q,uj m ;q',uj m/ )= 2 (2tt) 5 (q + q ) (3h5 m - m < 

n (iu m + hq 2 /2m)(-iu m : + hq' 2 /2m) , 
h [ul + e B {q) 2 /h ][u 2 m , + 6 B {q') 2 /h } 

With Eq. ( |3.51| ), we find the following general formula for the depletion in a disordered 
Bose superfluid, i.e., the reduction in the order parameter magnitude due to disorder, quan- 
tum and thermal fluctuations, 



d d q 
(2tt)° 



1 + Mqf 

U[q) ^ e 0e B (q) _ I 



n 1 ij^+^r (j^rn + W /2m)(-iu; m , + fig 2 /2m) 1 

71 fi P m y K + eB(g) 2 /^ ]K' + eB(<?) 2 A J J 

This expression depends on imaginary time, when the disorder breaks translation invariance. 
For a statistically translation-invariant system ( |4.7| ), the disorder contribution in Eq. ([4.1 8|) 

becomes 

n f d d q 1 ^ (hq 2 /2m) 2 - u 2 . , 

which for static defects ([4.1 3|) simplifies further to yield the temperature-independent result 

^^/t^^VaW. (4-20) 



4m 2 J (2ir) d \e B (q) 



These expressions will be evaluated explicitly for interesting cases in Sec. [IV C| . But we 
already see that for isotropic static disorder there is an intimate relation between the defect 
contributions to the depletion and the normal-fluid density, which are both temperature- 
independent, namely [ |TE| . [H| 

d 

mn A = -^PnA ■ (4.21) 

We conclude this section by deriving the effect of disorder on the normal-fluid density, 
using Landau's arguments. This analysis yields the most general expressions for the tilt mod- 
ulus tensor; note that Eq. ([4.15|) requires a statistically translation-invariant system (which 
is the most common situation), while the elegant derivation using affine transformations in 
App. [A] |H| is directly applicable only in the zero-temperature limit (/3h — > oo). 

For a superfluid in motion with respect to the surrounding walls and the intrinsic disorder, 
we have seen that the corresponding Green's function is simply obtained from the one for 
the system at rest by shifting the Matsubara frequencies according to u m — > oj m — iq ■ v 
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and uj m i — > ijj m ' — iq' ■ v [cf. Eq. (|3.63|) 1. Because the disorder correlator in the comoving 
reference frame reads 



A(q,a; m ;q',u; r) 



A(q, uj m + iq ■ v; q', u m > + iq' ■ v) , 



(4.22) 



these transformations imply that the disorder correlator to be inserted is effectively the one 
at rest. Thus, 
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A(q, u m ; q',^ 



(4.23) 



Upon expanding Eq. ( |4.23|) to second order in q - v and q' ■ v, and making use of Eqs. (|3.61 ) 
and ( 3.60D , we find for the normal-fluid density, which due to the disorder has become a 
tensor and in general a function of r as well, 
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(4.24) 



If translational invariance is statistically restored by averaging over the disorder, the defect 
contribution to Eq. ([4.24 ) reduces to the equivalent of Eq. (|4.15 ), because all the terms 
in the numerator of ( |4.24j ) which are linear in uj m then vanish. For an isotropic system, 
the normal-fluid density becomes a scalar, see Eq. ( |4.12| ). For possibly anisotropic static 
disorder, the renormalization of the normal fluid density becomes 



(5h 2 

-Ad 5 - 



d d q 



(27r) d Vsinh[/?e B (g)/2] 
which generalizes Eq. ( |4.14 ). 
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no 



d d q 
m J (2tt) c 



qiqj 



' h 2 q 



A(q) , (4.25) 



C. Results for point like, linear, and planar defects 

We now evaluate the density correlation function and the effective contributions to the 
depletion and the normal-fluid density caused by pinning to static disorder explicitly for 
the cases of point and (spatially) extended defects. We work in the boson picture, and 
use the words "point" , "line" , and "plane" to indicate the dimension of the disorder in 
the c?-dimensional spacelike direction r. According to Eqs. ( |4.20| ) and ( |4.25| ), the defect 



contributions for static disorder are independent of boson temperature; i.e., in the case 
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of periodic boundary conditions there will be no finite-size corrections in the imaginary- 
time direction for the corresponding quantities in type-II superconductors. For a flux liquid, 
such "static" disorder translates to defects correlated along the magnetic-field direction, and 
"point" and "linear" defects in d = 2 map on columnar pins and planar disorder, respectively 
(see Table |T| for an explicit comparison in d — 2). We shall return to these CclSGS, clS well 
as to "true" (uncorrelated in r and r) point disorder and tilted extended defects in Sec. [VT|, 
where we shall also compare the results obtained for periodic and open boundary conditions. 
For boson point defects, the disorder correlator reads 

A(r, r; r', r') = A 5(r - r') , A(q, uo m ) = A f3M mfi . (4.26) 

This immediately yields the density correlation function 

n hq 2 /m ( n h 2 q 2 ^ " 



ul + e B (q) 2 /h \me B {q) 2 ) 



and hence for the static structure factor 



— — n h 2 q 2 (3e B (q) A f n Q h 2 q 2 \ 



see also Eq. 

Furthermore, 

d A n r d d q ( hq \ 4 . m 3 nn d ^ _1 a (d//2) ~ 2 r°° x d ~ l 

m fTX — — ~p^~a = A — — / - — ^-7 ( — 7-— ) ^A ~r \ — dx , (4.29) 

4^" A 4m J (2n) d \e B (q) J 2 5 - d n 2 T(d/2)h d Jo (1 + x 2 ) 2 ' 1 ; 

where in the final equation V(q) ~ Vq = Aira/m was assumed (x = hq/Ay/nano). In two 
dimensions, one finds 

d = 2 : mnX = p^X/2 = Am 3 /16n 2 h 2 a , (4.30) 

while in d = 3 |T7|,|T§ 

d = 3 : mnX = 3p^/4 = Am 3 n /2 /8n 3/2 h 3 a 1/2 . (4.31) 

We can readily generalize these results to the case of static parallel defects, extended 
along a du -dimensional subspace rn. The disorder is therefore uncorrelated only in the d± 
transverse spatial directions rj_, i.e., 

A(r,r;r / ,r , ) = A(J(rj.-rl) , A(q, oj m ) = A (27r) d|l <5(q||) (3h5 mfi . (4.32) 

The dimensions d\\ = 1 and thus d± = d — 1 correspond to line defects, while d\\ = 2, 
d± = d — 2 describe planar disorder, etc. The corresponding dynamic structure factor then 
reads 

n Q kq 2 /m ( n Q h 2 q 2 ± 



S( q ,u m )= 2 7""/;'; /ft 2 +^l "V^ l (2vr) d il5(q||)/?M m , , (4.33) 
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and the static structure factor becomes 



n h 2 q 2 (3e B (q) ( n h 2 ql 



S(q, t = 0) = -^f- coth + A ( '•»" ^ 1 (27r) ci n5(q | |) . (4.34) 



The disorder-induced depletion follows from Eq. ( |4.29|) by replacing c? with d 



n 



f hq ± V m 2 nS d±/2)_1 a^/ 2 )- 2 /•<» x^" 1 



nA A 4m2i (2tt)^ Ub(?±)J A 25-^x^(^/2)^ X (l + x 2 ) 2 ^' (435) 

The normal-fluid density, defined as the transport coefficient characterizing the linear re- 
sponse of the system to motion with respect to its boundaries, will now depend on the 
direction relative to the disorder. If this motion is enforced along the disorder direction rn, 
the defects will have no effect [see Eq. (|4.25|) with qi = q\\. 



IMZ = ; (4.36) 

on the other hand, if the probe acts transverse to the defects, we find a similar result to 
Eq. ( f4.29| ) (for point disorder, all directions are transverse to the defects), namely 

m«A = jSIa- (4.37) 
For linear defects in two dimensions or planar disorder in d = 3 (d± = 1) this yields 

d± = l : mnX = P^Ta/4 = Am 3 /647r 3/2 hn /2 a 3/2 , (4.38) 
whereas for line defects in three dimension 

d± = 2 : mn£ = p^Xa/2 = Am 3 /16n 2 h 2 a , (4.39) 



which is, of course, identical to the result ( |4.30|) for point defects in d = 2. 



In the framework to the Bogoliubov approximation, which amounts to an expansion with 
respect to both thermal and quantum fluctuations, the disorder-induced depletion for d± < 4 
is smaller than the corresponding enhancement of the normal-fluid density, Eq. (|4.37| ). The 
bosons will become localized when p s ± — > (i.e., the corresponding "tilt" modulus c\ 
diverges !). The emerging new low-temperature phase of localized bosons is the Bose glass 
p0^p7[ (we use this term in a generalized form which also includes localization by extended 
defects), and is characterized by a vanishing superfluid density p s = 0. In the corresponding 
"Bose glass" phase of localized magnetic flux lines, the tilt modulus transverse to the defect 
directions diverges, which implies the "transverse Meissner effect", namely the absence of 
any response of the superconductor towards transverse magnetic fields (up to some critical 
tilt angle), see Sec. [VI F| . 

As a final investigation of how disorder affects superfluids, consider a Helium film on 
a substrate with disorder in the form of lines (deposited, say, via microlithography) with 
random positions and orientations. Imagine that the lines are arrayed in "families" indexed 
by a two-dimensional unit vector n perpendicular to every number of a family. Within a 
single family, the Fourier-transformed disorder correlator must have the form 
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Aft(q, u m ) = 2A5«(n • q) (3h5 m>0 , (4.40) 

where the one-dimensional delta function restricts only the wavevector component parallel 
to n. Upon using the integral representation of the delta function, 8(x) = /f^ e lxs ds/2ii, we 
can compute the quenched average over families in different directions by integrating over 
n = (cos 0, sin 0) to find 

1 r 27r A 
A(q,u m ) = — / d(/)A, h (q,u} m ) = — /3h6 m>0 . (4.41) 
Z7T jo q 

Results for the density correlations, depletion and normal-fluid density can be obtained by 
inserting this correlator into Eqs. ( |4.8| ), ( [4.20|) , and (|4.14|) ; thus we find, by deviding the 
corresponding results (|4.33|) , ( |4.34|) for point defects by q, the density correlation function 

n hq 2 /m A ( n h 2 q± 
^m + e B{q) 2 /h 2 q \me B (q A 



S (l^rn) = - 9 ./.,,» /t2 + ~ [ ^;_ ( :\ 2 (27r) d ii5(q||)/3M m , , (4.42) 



and the static structure factor 



n h 2 q 2 (3e B (q) A ( n h 2 ql x " 



S(<hT = 0) = coth^^ + - u ^ (2vr) d n5(q||) . (4.43) 

2me B {q) 2 q \me B {q±) 2 J 

Similarly, we find in analogy with Eq. ( 4.29| ) 

_d no r d d q 1 / Tiq \ 4 m 3 n d ~ 3)/2 a^- 5 ^ 2 r°° x d ~ 2 

mn A --p nA -A — J — -^—j « A ——^——^ ^ ___ rf x , 



(4.44) 

which is logarithmically divergent in d — 1, but yields a finite result in two dimensions, 

d = 2 : mn£ = p^X/2 = Am 3 /128vr 3/2 n /2 /ia 3/2 , (4.45) 

whereas in d = 3 

d = 3 : mnX= 3p^X/4 = Am 3 /167r 3 ^ 2 a . (4.46) 

We see that, although the correlator (|4.40| ) is singular as q — » 0, the disorder contributions 
to Eqs. (|4. 14|) and ( |4.20|) remain infrared-convergent (for d > 1). The depletion and normal- 
fluid density are, however, increased due to the correlations, relative to the same amount of 
disorder scattered about in a "point-like" fashion. 



V. DISORDER CONTRIBUTIONS FOR FLUX LIQUIDS (OPEN BOUNDARY 

CONDITIONS) 

We now investigate the influence of disorder on the correlation functions, "boson order 
parameter" , and the tilt modulus for flux liquids with open boundary conditions along the 
magnetic-field direction, using the results obtained in App. [Cj In this and the following 
section, we employ the flux line language, but retain the boson superfluid notation, in order 
to facilitate comparison with the results for the case of periodic boundary conditions derived 
in Sec. |TV| We shall discuss the equations for general disorder correlators, which lead to 
rather lengthy and not very transparent formulas. Specific results for a number of interesting 
types of defects are tabulated in Sec. 0. 
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A. Density and tilt correlations 

Upon carrying out the matrix products in Eq. ( |C9|) , we arrive at the following explicit 
results for the correlation functions with open boundary conditions in harmonic approxima- 
tion, to first order in the disorder correlator, 



(0(q, u m )Q(q',uj m/ )) 
me B (q) 2 /noh 3 q 2 



l/h 2 



a*. x 2h + u m u m , + e B (q) 2 /fi 2 \ , 
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e e(g) 



K + e B (g) 2 /^][^ + eB(g') 2 Al 
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2 wa + ee^Vft 2 , 
A(q, u; n ; q', uv) , (5.1) 
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2 ^ + 6 B (g') 2 A 2 
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tanh 



x bm'n' + 
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Pe B (q>) 



tanh ■ 



uJ 2 n , + e B (q>) 2 /h 2 



2 u; 2 n + e B (q) 2 /h 2 ) 
A(q,u n ;q',u n ') , (5.2) 



(7r(q,a; m )7r(q / ) a; m ')) 
riohq 2 jm 



^ m + e B (q) 2 /h 2 , 
n^q 2 q' 2 / m 2 



(3M ri 
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e B (?) ~ 2 ^, + e B (g) 2 A 2 7 

/, 2 /3e B (g) w m u; n - e B (q) 2 /h 2 \ 
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2 wa + ee^Vn 2 
A(q, o; n ; q', o; n /) . (5.3) 



PtB(q') Um'Un' - e B (q') 2 /h 2s 
2 ^, + e B (g') 2 A 2 



x <W,n' + 



^B(g') 



tanh 



These results should be contrasted with Eqs. Q4.4D - Q4.6D , which are valid for periodic bound- 
ary conditions, and may be recovered from ( |5.1| )-( |5T3| ) by taking into account only those 
terms in the brackets containing the Kronecker symbols 6 m - m ', <5 m ,n, and <5 m ', n '- 
From Eq. ( |5.3| ), we can deduce the equal- "time" density correlation function, 

"° ftV tanh^M(2 7 r) d % + q') 



S(q,r;q',- 



2me B (q) 



nlh 2 q 2 q' 2 



1 



E 



A(q,cu„; q',uj n > 



x 



m 2 e B (q)e B (q>) (^h) 2 [u 2 + e B {q) 2 /h 2 ][u 2 n , + 6 B (g') 2 /^ 2 ] 

" e B (g) ( -i« nT _ cosh[((3h-2r)e B (q)/2h] \ _ . sinh[(/3h - 2r)e B (q)/2h] 
h V cosh[/?e B (g)/2] J ^ n cosh[/3e B (g)/2] 

sinh[(/%-2r)e B (g , )/2fr] 



x 



x 



e_M (^ nlT _ cosh[(/m - 2r)e B (q')/2h}' 



cosh[/5e B (g')/2] 



cosh[/3e B (gO/2] 



(5.4) 
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generalizing both Eqs. (|4.9| ) and (|3.77|) . The ^-dimensional structure factor is now a function 
of r, implying that the mean-square variations of the vortex positions depend on the depth 
of the considered cross-section perpendicular to the magnetic field in the sample. On the 
surfaces at r = and r = j3h, the fluctuations vanish. Similarly, for the leading term of the 
current correlation function (i.e., the one originating in the phase fluctuations) we find 



C^(q,r;q',r) 



n m^Mtanh 



faB(g) p L,_ . ,.-,_„/, 

2 ij 



PtM) (27r) d 5(q+q') 



+n q t q j 



E 



A(q, u) n ; q',u; n > 



x 



ILO, 



W tt> K + ^(qymiul, + e B (q>) 2 /h 2 ] 

coshpft - 2r)e B (q)/2h]\ e B (q) sinhpfi - 2r)e B (q)/2h] 



cosh[/3e B (g)/2] J % 



cosh[/3e B (g)/2] 



x 



X 



, T _ cosh[(^-2r)6 B (g / )/2^] \ _ e B {q') sinh[(^ - 2r)e B (g')/2ft] 



cosh[/5e B (g')/2] J 



cosh[/5e B (g')/2] 



. (5.5) 



Eq. ( 3-30|) then yields the tilt correlations. 

In the center of the sample, r = /3H/2, these expressions simplify considerably; upon 
assuming translational invariance in the d transverse "spatial" directions, the density and 
current correlators read 



S(q,r = /?ft/2) = tanh^ 4 



x e 



2me B (g) 
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1 
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-iu) n f3h/2 



(5.6) 



1 



cosh[/3e B (g)/2] 
ragg 2 A(q, u; w ; -q, ay) a> n a; n / 
( ;//)^O^K + e B ( ? )V^][^, + e B (g)V?i 2 ] 
1 



E 



cosh[/3e B (g)/2], 



cosh[/3e B (g)/2] < 



itf(q). (5.7) 



Note that e _ For static defects ( |4.13j ), the disorder contribution to Eq. ( p77[ ) 

vanishes, while the static structure factor becomes 



— — — - n h 2 q 2 /3e B (g) 

5(q, r = = ^"^y tanh _ + A(q) 



2^2 



me B (qf 



1 



1 



cosh[/3e B (g)/2] 



(5.8) 



B. Depletion and vortex tilt modulus 

Within the Gaussian approximation (|3.42|) , (|5 . 1| ) — ( |573D we can readily obtain the disorder 
contribution to the Green's function for open boundary conditions, supplementing the pure 
part ( gjp , 
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'ftq 2 



n /ft 2 



K + eB(q) 2 /h 2 ]iul, + e B ( q >y/h 2 }t;, 
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+ iu m b m ^ n + 
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u 2 n , + e B (q>) 2 /ft 2 



(5.9) 



Upon assuming as usual translational invariance in r, this leads to the disorder-induced 
depletion, i.e., a disentanglement of pinned vortex lines, 
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d d q 
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(5.10) 



compare Eq. (|4.18 ) for periodic boundary conditions. It is a straightforward, but rather 
tedious task to compute the disorder renormalization of the normal-fluid density from 
Eqs. ( |3.60|) and ( |3.61 ), following Landau's arguments (see Sec. [IV B|) ; besides expanding 
to second order in q • v, this involves a number of Matsubara frequency sums for which the 
list in App. |D 1| is helpful. The final result reads 
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(5.11) 



and is to be contrasted with Eq. fl4.24| ). Note that as a consequence of the open boundary 
conditions (which suppress fluctuations in the boson order parameter at the boundaries), 
both results ( |5.10| ) and ( |5.11| ) vanish for r = and r = /3H, respectively. At r = /3ft/2, these 
expressions simplify, and the depletion and normal-fluid density in the center of the bulk 
become 
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(5.12) 
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For static disorder (|4.13|) , we find the following comparatively simple expressions 

^a(t" = pft/ 2 ) = 
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(5.14) 



., A (,^ft/2)^/|^^(-^) 2 A( q )x 



43 



x (i 1 ^(g)/2 /?e B (g) \ 

^ cosh[/?e B (g)/2] sinh[/3e B (g)/2]^smh[^B(g)]y ' l ' J 



In contrast to Eqs. ( |4.2Q|) and ( |4.25| ), these results do depend on (3%. With the defects being 



"oriented" along the r direction, it is plausible that there are no finite-thickness correc- 
tions to the disorder contributions for quantities such as the depletion and the normal-fluid 
density, provided periodic boundary conditions are applied; for all points along the de- 
fect (world) line are then fully equivalent. This explains in physical terms why the effects of 



static disorder on Bose superfluids turned out to be independent on temperature, see Sec. [TV. 
With open boundary conditions, on the other hand, the enhanced thermal wandering near 
the surfaces comes into play (compare Sec. [11 B|), which renders the defect influence depen- 



dent on the sample thickness. As is apparent from Eqs. Q5.6|) — fl5.8|) and (|5.12|) — (|5.15|) , these 
finite-size contributions typically decay exponentially with growing sample thickness /3h, the 
characteristic length (the analog of the thermal de Broglie-wavelength for superfluids) being 
cx /3hci at small wavenumbers. Note in addition that there appears no simple relation such 
as ( |4.21| ) or (|4.37|) between the disorder contributions to the depletion and the normal-fluid 
density at any (3% < oo any more. We defer the explicit evaluation of Eqs. (|5.12|) and ( |5.13| ) 
for several types of defects to the following section. 



VI. EXPLICIT RESULTS FOR VARIOUS KINDS OF DISORDER 

We now specialize the cumbersome expressions obtained for general disorder in the previ- 
ous section to specifically interesting types of defects, and in some cases list the corresponding 
results for both periodic and open boundary conditions along the r ~ z direction. Although 
we have the application to magnetic flux liquids in high-temperature superconductors in 
mind, we express our results in the boson superfluid notation in d "transverse" dimensions. 



Only in the final subsection [VI F| , when we summarize a number of experimentally relevant 
consequences of our investigations, we put d = 2 explicitly and translate to the flux line 
language. 



A. Uncorrelated disorder — point defects 

We start with disorder completely uncorrelated in r and r, 

A(q, uo m - q', O = A (27r) d 5(q + q') /?M m ,-m' ; (6.1) 

the quenched disorder average restores translational invariance in "space" and "time". In 
the flux line picture, Eq. ( |6.1| ) describes point defects, e.g. oxygen vacancies in the cuprate 
planes. (For "true" bosons, a realization of this space- and time-dependent disorder is 
more difficult to conceive.) For this highly "dynamical" disorder ( |6.1| ), we find the following 
defect contributions to the structure factor and the leading component of the "static" current 
correlation function (stemming from the phase fluctuations only), first for periodic boundary 
conditions [see Eqs. ( P~9|) and flQpl, 
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while the results for for open boundary conditions [Eqs. (|5] 
the sample (r = /3h/2), read 
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As in the pure case, the limit f3Ti — > oo is thus approached from different sides for periodic 
and open boundary conditions, respectively. 

In order to compute the depletion from Eqs. Q4.19Q and (|5.12|) for periodic and open 
boundary conditions, we use the Matsubara frequency sums tabulated in App. p 1| . The 
results are 
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(6.7) 



In the limit /3ft — > oo, these expressions coincide, of course, 
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With x = hq/4y/irnoa, we find for a delta-function repulsive potential [V(q) ~ Vq = Ana/m] 
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in two dimensions, this yields 
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and in the three-dimensional case 
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By exploiting the identities f£°[x d 2 (cothx — 1 
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1)] J™ x 4 - 1 / sinh 2 xdx and / °°[x*- 2 (l - tanhx) - x d - 1 / cosh 2 x]dx = -[(d - 2)/(d - 1)] 
/ °° x d ~ l / cosh 2 xdx for d > 2, we can derive the leading finite-temperature / sample thickness 
corrections in the phonon approximation, 63(9) ~ ^cig (x = (3hciq/2), 
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with Eqs. (P15|) and ( P16| ) of App. P 2| we finally find in two dimensions 
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e.g., in three dimensions, 
d = 3 : 
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(6.14) 
(6.15) 

(6.16) 
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(6.18) 
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As expected, the sign of the finite-temperature / size corrections is opposite for periodic and 
open boundary conditions, respectively, again reflecting the fact that the enhanced thermal 
wandering near the surfaces increases the tendency for entanglement of the (world) lines in 
the case of open boundary conditions. In d = 2 the leading correction oc k B T vanishes for 
open boundary conditions. 

With somewhat more effort, the disorder renormalization of the normal-fluid density 
/ tilt modulus may be obtained evaluating ( |4.12|) and (|5.13|) , respectively; as the disorder 
correlator is isotropic, one finds 
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n /? 2 ft 4 f d d q Aq 4 cosh[/3e B (g)/2] 
8dm J (27r) d %e B (q) sinh 3 [/3e B (g)/2] 



p.b.c. : p nA = ^3— / ^ —-—3— - — - , (6.20) 



, n ^h 4 f d d q Aq 4 1 

° ' P " A " 8dm J (27r) d ^ B (g)sinh[/3e B (g)] ' 1 } 

note that these results for periodic and open boundary conditions are not simply related, 
because the normal-fluid density arises from the next-to-leading order in a long-wavelength 
expansion. However, in the limit (3Ti — > oo one has 



Pua(T = 0) = , (6.22) 

implying that mere point disorder cannot alter the tilt modulus of a flux liquid in thick 
samples; see also App. |A|. For finite "temperatures", the leading corrections become for 
periodic and open boundary conditions, respectively [using J oo x d+2 (coshx/sinh 3 a;)(ix = 
\{d + 2)/2] Jo 00 x d+l / sinh 2 x and Eq. flDT5p 1 

nbc • a — m T (d + S)C(d+l) An (k B T)^ 

p.b.c. . p„ A (T)- 2d+l7rd/2T{1 - d/2) > ( 6 - 23 ) 

1 \ T(d + 3)C(rf + 3) An (fc B T) d+1 



o.b.c. : Pn^T)=-(l-^ 2d+W{i + d/2) -^-p , (6.24) 
which in two dimensions becomes 



d = 2: p.b.c: p nA (T)= 3C(3)An (fc B T) 3 /7rm^ , (6.25) 
o.b.c. : p nA (T)= -93C(5)An (k B T) 3 /6ATrmh 3 cl , (6.26) 



while in three dimensions 



d = 3: p.b.c: p nA (T) = 7v 2 An (k B T) 4 /9mh 4 cl , (6.27) 
o.b.c. : p nA (T) = -7r 4 An (A; B T) 4 /192m^ 4 c[ . (6.28) 

Again, these finite-temperature / sample thickness corrections display opposite signs for the 
two different boundary conditions. 



B. Correlated random disorder — nearly isotropic splay 



For magnetic flux lines in high-temperature superconductors, a deliberate introduction 



of splayed columnar defects has been suggested as a means to enhance flux pinning [44|,f45 



In the flux line notation, the disorder correlator for such randomly tilted columnar defects, 
each described by a trajectory Ti(z) = R; + ViZ, with a Gaussian distribution of the tilts Vj, 



-P[ v i] n« e Vi ' 2vD , reads in Fourier space H,fR: 



A(q,g 2 ) 



Ai 



-q 2 j2v 2 D q 2 



2nvoq 



(6.29) 



Upon taking the limit Vd — > oo, with Ai/vp held fixed, one has A(q,q z ) oc 1/q, i.e., an 
"almost" isotropic situation; as this correlator was obtained by starting from a distribution 
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centered around the z axis, which thus constitutes a preferred direction, this limit does 
not quite lead to the "truly" isotropic limit A(q, q z ) oc l/(q 2 + q 2 )" 1 ^ 2 [compare the two- 
dimensional correlator ( fj.41|) 1- The additional q z dependence here is not expected to affect 
the physical implications drastically, but would render the evaluation of the depletion and 
tilt modulus considerably more cumbersome. 

We thus study the simpler case of "neaarly isotropic" splay, namely 



A(q,u m ; q', u. 



ml 



A 

q 



]2ir) d 5(q + q') /3h5„ 



(6.30) 



With this correlator, we can immediately take over the results for point defects from the 
previous subsection, because none of the Matsubara frequency sums is changed; we simply 
have to divide A by q everywhere. Hence 



p.b.c. : SA(q) 



A 



2„2 



n Q h q 



% e B(<?) \2me B (q) 



coth 



CijA(q) 



A 



n$Tiq 



%£B(g) V 2 
for periodic boundary conditions, while 

;2„2 \2 



coth 



Pe B (q) 
2 

Pe B (q) 



1 + 



1 - 



o.b.c. : Sa(<i)-- 



A 



CyA(q) = 



' n h q 2 
%e B (g) V2me B (g) / 



n Q hq s 



tanh 



Pe B (q) 



A 



hqe B (q) 



smh[(3e B (q)} 

Pe B (g) \ 
sinh[/3e B (g)]/ 

Pe B (g) 
sinh[/3e B (g)] 

Pe B (g) \ 
sinh[/5e B (g)]/ 



(6.31) 
(6.32) 



^(q) 



(6.33) 
(6.34) 



for open boundary conditions, evaluated at r = (3%/ '2. 

In the same manner, dividing the integrands of Eqs. ( p.6| ) and (|6.7|) by q yields the defect 
contribution to the depletion 



p.b.c. 



no 
4 



d d q 

(2^F hqe B {q) 

2^2 \ 2 



A ooth^Mx 



X 



h z q 2 



Pe B (q) 



Pe B (q) 



o.b.c. 



n A -- 



no 
4 



d d q 



2me B {q) 
A 



sinh[/5e B (g)] 



(27r) d Hqe B (q) 

2^2 \ 2 



tanh^Mx 



x 



h z q 2 



(3e B {q) 



2me B (q) 



smh[(3e B (q)] 



In the limit (3% — > oo, both expressions lead to 



n 



d d q 



A 



4 J (2ir) d hqe B (q) 



2„2 



h 2 q 



2me B (q) 



-A 



- 1 



sinh[/5e B (g)] 



Pe B (q) ' 
sinh[/5e B (g)] 



d d q qV(q) 



(6.35) 



(6.36) 



4m J (2vr) d e B (g) 3 ' 



(6.37) 



For V(q) ~ Vq = 4ira/m (x = hq/ 'A^nano) we thus find for the disorder contribution to 

n - = n - |(^)| 2 , 
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-A 



mn 



(d-l)/2 _(d-3)/2 



3-d 



d-3 



7T 3 / 2 r((i/2)^ 7 (1+X 2 ) 3 / 2 



■ dx 



(6.38) 



which is negative and logarithmically divergent in two dimensions, suggesting that KV')! 2 
becomes unbounded above! Evidently the flux lines are easily trapped and swept along by the 
disorder; thermal wandering presumably becomes superdiffusive, leading to the divergence of 
TTX (see Ref. |32| for a related case). The mapping onto nonrelativistic bosons is inadequate 
to describe completely this "sup erent angled" state of interacting lines. In three dimensions, 
the corrections are finite, 



2fc3 



n A = —Amno/An % 



(6.39) 



and the wandering remains diffusive. Upon dividing the expressions of the previous subsec- 
tion by q, we find for the finite-temperature / sample thickness corrections for d > 2, 



p.b.c. : 
An A (T) 

o.b.c. : 

An A (T) 



An (k B T) 



d-2 



87r d / 2 r(d/2)^cf- 1 
An (k B T) d - 2 



Sd,3 + 



d-3 f°° x d ~ 2 



d — 2 jo sinh 2 x 



dx + 



{k B Tf 



x 



m 2 Ci Jo sinh x 



-dx 



8n d / 2 T(d/2)h d c d l - 1 
which for d — 3 yields 



d-3 



x 



d — 2 jo cosh x 



dx + 



(fc B T) ; 
9 4 



.r 



cosh x 



dx 



d 



p.b.c. : A^(T) = + 0(T 3 ) 



o.b.c. 



Aix 2 h 3 cl 
Atz a (T) = -0(T 3 ) 



while for d > 3 



p.b.c. : An A (T) = 



(d - 3)r(d - 2)C(d - 2) An (fc B T) 



d-2 



o.b.c. : An A (T> 



2<W 2 r(d/2) ^cf 1 

1 3)r(d - 2)C(d - 2) An (k B T) d ~ 2 



2 d-3 



d „d-l 



In the same manner, we obtain for the normal-fluid density, 

n (3 2 ^ f d d q Aq 3 cosh[/3e B (q) /2] 

p.b.c. : p„ A - ' 

o.b.c. : p^X~- 



8dm J (2ir) d he B (q)smh 3 [f3e B (q)/2] 
n (3 2 h 4 f d d q Aq 3 1 



8dm J (2ir) d he B (q) smh[(3e B (q)] 



(6.40) 
(6.41) 

(6.42) 
(6.43) 

(6.44) 
(6.45) 

(6.46) 
(6.47) 



which approaches zero as j3h — > 00, as is shown to be the case for any disorder correlator 
A(q, u m ) independent of oj m in App. [A]. For finite temperatures or sample thickness, 

T(d + 2)C(d) An (k B T) d 



p.b.c. 
o.b.c. 



PnA(T)-- 



PnA(T)-- 



2 d + i 7r d/ 2r ( 1 + d / 2 ) m n d c d+3 ' 

1 \ T(d + 2)C(d + 2) An (k B T) d 



1 



2 d+2j 2 d+2 7[ d/2 T ( 1 + d / 2 ) mh d c d 



d d+3 ' 



(6.48) 
(6.49) 
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in two dimensions this becomes 



P-b.c. : p nA {T)= 7rAn {k B T) 2 /8mh 2 cl , (6.50) 
o.b.c. : p„ A (T)= -n 3 An (k B T) 2 /256mh 2 cl , (6.51) 



and for d = 3 



d = 3: p.b.c. : p nA (T)= 2C(3)An (fc B T) 3 /v 2 m/i 3 c? , (6.52) 
o.b.c. : p nA (T)= -31C(5)An (A;BT) 3 /327r 2 m^ . (6.53) 

We see that on/y for d = 3 and with open boundary conditions there exists a simple relation 
between the finite-size corrections to the disorder renormalization of the depletion and the 
normal-fluid density, namely 



d = 3 , o.b.c. : p nA (T) = [31C(5)/9C(3)] mAn A (T) . (6.54) 

Finally, we remark that in the above calculations we have not taken into account the 
possibility of a direct coupling of the defects to the tangent field (i.e., the momentum or 
current density in the boson language), but have exclusively studied the effects of disorder 
which couples to the particle density, i.e., changes the chemical potential locally. The latter 
mechanism should remain the most important effect of the disorder even in the case of 
splayed correlated defects; however, it is conceivable that in addition a term of the form 
t£)(r, r)-g(r, r), see Eq. (|2.49|) , with t£»(r, r) representing another quenched random variable, 



may have to be included at least in an "effective" theory for which one further step of some 
course-graining procedure has been carried out. We neglect this scenario here, as well as in 
Sees. |VTDl and [VTEl below. 



C. Correlated disorder — parallel untilted extended defects 

We now turn to the analog of static disorder for bosons, namely untilted extended defects 
in flux liquids. We can use the correlator ( |4.32 ) , with the understanding that dii = 



describes point defects in the boson representation, i.e., linear (columnar) defects in the flux 
line picture, and similarly d\\ = 1 represents defect lines / planes, respectively, etc. The 
results may then simply be expressed in terms of the transverse dimensionality d± = d — d\\. 

The disorder contribution to the static structure factor with periodic [see Eq. (|4.34j )1 and 
open boundary conditions [cf. ( |5.8|) 1 read 

n h 2 q]_ 



p.b.c. : fif A (q)= A u ^ (2vr) d n ( 5(q||) , (6.55) 

Vme B (g±)V 



n h g ± 



2„2 \ 2 / i \ 2 



whereas the leading part of the current (tilt) correlators vanish in this case in the static 
limit. 

The results for the depletion and the normal-fluid density in the limit (3h — > oo have 
already been explicitly given in Sees. |IVB| and [IV C| , where we had also seen that for periodic 



boundary conditions the finite-temperature / sample thickness corrections vanish entirely, 
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p.b.c. : An A (T) = 0, p n ijA{T) = . 



(6.57) 



However, for open boundary conditions there are finite-size corrections, induced by the 
stronger thermal fluctuations in the vicinity of the sample surfaces. Eq. ( |5.14j ) yields (in the 
sample center, r = (3%/ 2) 



o.b.c: An l(I1 = -A^/ 



n r d d± q± ( hq^ 



1 



2m 2 J (2vr)^ \e B (q ± ) J \cosh[(3e B (q) /2] 2 cosh 2 [/3e B (g)/2] , 



(6.58) 



which can be evaluated in the phonon approximation using Eqs. QD18 ) and ( P16j ) of App. |D 
where also the function £ (d) is defined, 



o.b.c. : Aua(T) 



2T(d 



± 



TT d ^ 2 T(d ± /2) 



ad±) - z- d± (l - 2 2 -^) C(d± - 1; 



An (k B T) d ^ 

(6.59) 



for the relevant transverse dimensionalities this reduces to [£(2) = G ~ 0.916 denotes 
Catalan's constant] 



<i_|_ 
<i i = 3 



1 : An A (T) = -(tt - l)An Q k B T/2Tmi 2 hc\ 
2 



An A (T) = -(AG - ln2)An (fc B T)727rmVc? 
An A (T) = — (37T - l)A7i (fc B T) 3 /12m 2 7i 3 c£ . 



(6.60) 
(6.61) 
(6.62) 



The normal-fluid density along the direction of the extended defects remains unaffected 
by the disorder, 



o.b.c. : p n || A CO = . 
In the perpendicular directions, Eq. ( |5.15| ) leads to 

An f d d± q± ( hq± \ 4 



(6.63) 



o.b.c. : p„_laCO = 



d ± mJ (2vr)^ \e B (q ± ) 
( 1 



x 



(3e B (q A 



cosh[(3e B (q ± )/2] sinh[/3e B (g±)/2] sinh[/3e B (g±)] 
which in the phonon approximation is readily evaluated using Eqs. (P17 ) and (|D18|) 

2T(d ± ) 



, (6.64) 



Pn±A(0 : 



vr^/ 2 r(i + d±/2) 



1 - 2- d A (1 - 2- ( ^ +1) ) d ± ((d ± + 1 



An (k B T) d± 

mh d ^c d±+i 



(6.65) 



Explicitly one gets 
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d± = 1 : Pu±a(T)= -(tt + 4)An fc B T/4m^ 

= [7r(7r + 4)/2(tt - 1)] mAn A (T) , (6.66) 
d± = 2 : p n±A (T)= - [16G + 21C(3)] An (k B T) 2 / 'Mmh 2 c\ 

= [(16G + 21C(3)) /(16G - 4 In 2)] mAn A (T) , (6.67) 
di = 3 : p„±a(T)= -7r(77r + 8) An (k B T) 3 / A8mh 3 c\ 

= [n(7n + 8)/4(3vr - 1)] mAn A (T) ; (6.68) 

thus the finite-temperature / size corrections to the disorder renormalization of the deple- 
tion and the normal-fluid density are proportional to each other for these parallel untilted 
extended defects, albeit with a complicated numerical factor. 



D. Correlated disorder — parallel tilted extended defects 

More generally, we can also discuss "moving" or tilted parallel extended defects, say, 
along the x direction; the corresponding disorder correlator then becomes 

A(r, r; r', r')= A 5(r ± - r' ± + v(t - r')e x ) , 

A(q,u m )= A (27rp5( qil ) /3h5 Wm , qxV , (6.69) 

generalizing Eq. ( (4.32Q . Recall that r_i_ refers to directions perpendicular to a set of tri- 
dimensional defects; (dn,d±) = (l,d— 1) means line defects, while (d\\,d±) = (2,d — 2) 
corresponds to planar disorder. For magnetic flux liquids, the tilt parameter is simply the 
tangent of the tilt angle, v = tana, see also Eq. ( |A12| ) in App. |A|. We shall restrict ourselves 
for simplicity to thick samples in this and the following subsection. If one wanted to apply 
periodic boundary conditions, one would have to keep the transverse momenta discrete as 
well, q x ^k = 2irk/L± and ensure that f3hv/L±_ = m/k is a rational number. 
In the limit /3H — > oo, Eqs. (|4.9|) and (|5.6|) yield the static structure factor 



2me B {q) \qp 2 + e B {q ± ) 2 /% J 

while the leading contribution to the static current correlation function [see, e.g., Eq. fl5.7|)1 

reads 



Cy(q) 



n me B (q) A ( n q ± q x v \ f n^\d [lX(n 
2 \q 2 v 2 + e B {q±y h J 



pIM) ■ ( 6 - 71 ) 



For v — *■ 0, these results of course reduce to those of the previous subsection. In real time, 
the density and phase correlation functions display singularities when v — > c±, i.e., when the 
defects move at the speed of sound of the superfluid. 

The disorder contribution to the depletion ( |4 . 1 9| ) becomes 

a ^0 f d d ^q± (hqj/2m) 2 - q 2 x v 2 

UA ^h 2 J (27r)^ [e B (q ± ) 2 /h 2 + q 2 x v 2 } 2 ' ^ ] 
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this can be easily seen to diverge for any nonzero tilt v in the cases d± = 1 and d± = 2, 
which we again interpret as indicating that the world lines are "dragged" along by the tilted 
extended defects, leading to super diffusive behavior. Only for v = does the depletion 
remain finite, see Eq. (|6.59| ). For d± = 3 (point disorder in a three-dimensional superfluid 
moving at constant speed v) one finds for delta-function interaction with strength Vo = 
4na/m, 



d, 



Am n. 



2 1/2 



o 



1 



(6.73) 



8tt3/^V/ 2 (l + mV/47m a) 1/2 ' 

The disorder renormalization of the normal-fluid density and tilt modulus perpendicular 
to the extended defects can be inferred from Eq. fl4.15| ), see also App. [A], 

jd\ „ ^ f„ \2 /fc2 o^,2„,2 



Pn ±i A 



A 



no 



d d ^q± 2 2 e B (q ± ) 2 /h - 3q 2 x v 
Q±Q± 



(2tt 



while p n \\A = 0, of course. In the case d± ■ 
and 

Am 3 



P, 



1, there is only one perpendicular direction, 



m v /2-Krioa 



i . - • (6 . 75) 

167r 3 / 2 /m a 3 / 2 (1 + m 2 v 2 /A-irnoa) 

for d± = 2, on the other hand, we have to distinguish between the direction along the tilt 
and perpendicular to it, 

1 + m 2 v 2 /27rn a 



d. 



Pn±xA z 



Pn ±yA z 



%h 2 v 2 

Amno 
irh 2 v 2 



[1 + m 2 v 2 1 'Aim^a) 
1 



3/2 



1 



1 



and finally, for d 
d i = 3 : 



Pn±xA z 



Pn±yA z 



3 one gets 
4An 2 1 a 



h 3 v 3 

2An 2 a 

h 3 v 3 



arsinh 



mv 



m?v 2 /Airnoa) 1 ^ 2 
mv / (Annoa) 1 / 2 



(47rn a) 1 /2 ^ + m 2 y 2 /4 7rnoa ) 1 /2 _ 



(47rn a) 1 / 2 



1 + 



2 2 \ 1/2 



arsinh 



mv 



(47rn a) 1 / 2 



(6.76) 
(6.77) 

(6.78) 
(6.79) 



These results show that for general nonzero tilts quite complicated formulas arise; and it is 
only in the very special situation v = that simple relations of the form p n ±A oc nX emerge. 

It is interesting to compute the fraction of particles / lines that are actually "dragged" 
along by the defects, 



1 



mvn v 



With Eq. (gTID we find 



n, 



A 



n 



phtt 

d d± q± 



(q = 0,^ m = 0)) 



q±q x 



(6.80) 



(6.81) 



m"J (2n) d ± \e B {q ± ) 2 /h 2 + q 2 v 2 J 

For v — ► this precisely coincides with the disorder-induced normal-fluid density 
[Eqs. ( |4.35|) , (fOTD l; for static disorder, the superfluid density is thus reduced by exactly 
the average fraction of particles that are pinned to the defects. 
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E. Two families of symmetrically tilted extended defects 



As a final step of generalization, we now introduce two families of extended defects, which 



are symmetrically "moving" (tilted) along the x direction The disorder correlator in 
this situation reads 



A(r,r; r',r'1 



A f 5(r ± - r' ± + v(t - r')e x ) + 5(r ± - r' ± - v(r - r')e x ) 



4 



+S(r ± - r' ± + v(t + r')e x ) + 8(r x - r' ± - v(r + r')e x ) 



A(q, u m - q', u m >)= — (2vr)^(q + q') (27r) d H<5(q||) x 

x (/%) 2 (5 m _ m / + 5 m , OT /) ((L™,^ + <L m • (6.82) 

It is important to realize that translational invariance along the r direction is broken here, 
which originates in the fact that there is a strong correlation between defects "moving" in 
the positive and negative x directions; or, equivalently, of defects "moving" forward and 
backward in imaginary time. As a consequence, we have to resort to our most general 
formulas in Sees. [TV| and [V], and Eq. (|A11|) , for example, is not applicable; but again, we 



shall only discuss the limit f3h — > oo here. Note that in the limit v — > 0, the results of 
Sec. [VI Q for untilted parallel defects are recovered. 

Interestingly, with this correlator the static structure function, see Eq. (|5.4|), for example, 



becomes identical with the result ( |6.70| ) for a single family of tilted extended disorder. This 
must be the case, as in a fixed cross section at r, those two cases look the same. However, in 
the leading contribution to the current correlation function, the effects of the two oppositely 
"moving" defect families precisely cancel each other, leaving behind just the pure result (this 
is actually true for finite (5h as well, with any boundary condition applied). This implies 
that in the Green's function, from which the depletion and the normal-fluid density are to 
be determined, we shall face a partial cancellation of the disorder influence, namely in those 
contributions stemming from the phase correlations, respectively. 

Again because of cancellations, the disorder-induced depletion, obtained from Eq. ( p,12| ), 



turns out to contain only the first term in the numerator of Eq. ( |6.72|) , 



— _ a _^o_ f d d ±q± ( q\ \ 

UA ~ A 4m2 J (2tt)^ \e B (q ± y/h 2 + ) ' 1 ' 

i.e., the leading contribution (in a long-wavelength expansion), which had caused the di- 
vergences discussed in the previous subsection, actually has disappeared. Therefore the 
depletion remains finite even for nonzero tilt. The explicit results read 

Am 2 1 

d± = 1 : n A = T-p. rrr , (6.84) 

647r 3 / 2 /mo /2 a 3 / 2 (1 + m 2 v 2 / 4rm af 2 

Am 2 1 

d± = 2 : ?2a= Ta7 i (6.85) 

16rfa (l +m 2 v 2/ 47rnoa )V2 ' ^ > 

d± = 3 : n A = 3 arsmh — — . (6.86) 

AnTfv (47m a) 1 / 2 



54 



In a similar manner, the disorder renormalization of the normal-fluid density is reduced, 
and Eq. ( |5.13| ) leads to 



pn±iA = A— / , u q±q ±i 
m J (2ir) a± 



e B (q±) 2 /h 2 - q lv 2 

[e B (g±) 2 A 2 + <7> 2 ] 3 



(27 



to be contrasted with Eq. ( |6.74|) . For d± — 1 this becomes 

Am 3 1 — m 2 v 2 /87rnoa 



d± = 1 : Pn_LA 



167r 3 / 2 ^no^ 2 a 3 / 2 (1 + m 2 v 2 /4irn a) 5 ^ 2 



(6.87) 



1 — m v /87rnoa 
1 + m 2 v 2 /Annua 

(6.88) 



For g?j_ = 2 we find 

rf± = 2 : Pn±x A = 



1 



Am 3 

37r 2 ft 2 a (1 + m V /47rn a) 3/2 1 + m 2 v 2 /A7rn a 
Am 3 1 



Pn±yA z 



M 2 h 2 a (l + m 2 v 2 /4nn a) 1/2 
and finally in three "transverse" dimensions, 

Anla 



2mn/ 



(6.89) 
(6.90) 



d± — 3 : p n ±xA 



Aula 



Pn±yA 



Ah 3 v 3 



h 3 v 3 

mv 



arsinh 



mv 



mv 1 — m 2 v 2 / Aim^a 



(47rn a) 1 /2 (4 7 m a) 1 /2 (i + ^2/47^0) 1/2 



(47rn a) 1 / 2 



1 + 



2 2 \ V2 



2 2 \ 
m v \ 



2iin Q a 



arsinh 



mv 



(47rn a) 1 / 2 



(6.91) 
(6.92) 



F. Experimental consequences for vortex physics 

We conclude by translating the results of the previous subsections into the flux line 
language, specializing to d = 2 dimensions, and discussing some of the experimental con- 
sequences of our results. As noted in Sec. fl B , we approximate the inter-vortex repulsion 



by V(q) = Vq/(1 + X 2 q 2 ), with V Q = (fil/An and the magnetic flux quantum (fi = hc/2e. 
As shown in Table [| the particle mass is to be replaced by the effective line tension 
€\ = (mj_/m z )eo ln(A/£), where the energy scale is set by eo = (0o/47rA) 2 ; we shall also 
need the scattering length 

a = e 1 K)/4vr = e 1 (0o/47r) 2 , (6.93) 

and the characteristic velocity 

ci = (n V?i) 1/2 = (n o 2 /47re 1 ) 1/2 . (6.94) 

With these substitutions, and the replacements h —>■ k^T and (3h — ► L, the structure factor 
and other correlation functions are readily obtained from our results in the previous sections. 
In the following, we shall review the explicit consequences for the fraction of entangled lines 
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as well as the renormalization of the tilt modulus by disorder and through modification of 
the boundary conditions. 

For a pure system, we found in Sec. [Tl that the mean "boson" order parameter squared 
n = \ (ip)\ 2 as function of temperature in the limit of thick samples L — > oo is 



n (T) 



n 



1 + a/(k B Tf 



n 



-i 



(6.95) 



(ip) measures the degree of entanglement, which increases with increasing temperature. The 
vortex contribution to the tilt modulus [see Eq. (|TTT|)1 is c^ -1 = (noei)" 1 . The leading 
finite-size correction for sample thicknesses L < oo to the fraction of disentangled lines 
diverges logarithmically, if periodic boundary conditions are employed along the magnetic- 
field direction, implying that in fact the "boson" order parameter vanishes in this case. On 
the other hand, for open boundary conditions, which are more realistic for most samples, 
the reduction in the "boson" order parameter is finite, cf. Eq. ( |3.80| ), 



An(L) = -(ln2)ei/27rfc B TL 



(6.96) 



To leading order in the sample thickness corrections, Eqs. ( |3.67|) and ( |3.86| ) yield for the tilt 
modulus with periodic and open boundary conditions, respectively, 



c\f l {L) = {nhY 1 1 =F 3((3)k B T/2nn~e 1 c 4 1 L 3 



p.b.c. 
o.b.c. 



(6.97) 



with C(3) ~ 1.202. As already discussed, the different signs for the finite-size corrections 
reflect the fact that for open boundary conditions thermal flux line wandering is effectively 
enhanced near the sample surfaces, and therefore both entanglement and response to exter- 
nal tilt are facilitated; the firm restrictions of periodic boundary conditions, on the other 
hand, reduce the number of entangled lines, and stiffen the system towards external tilt- 
ing. According to Eq. (|1.2|), any renormalization of the tilt modulus may also be viewed as 



a changing the mass anisotropy m±/m z < 1, with an increasing implying a smaller 

effective mass ratio and vice versa. Any enhancement of the normal-fluid density, and corre- 
spondingly a reduction of c\^~ x , either through finite-size effects as in Eq. ( |6.971 ) for periodic 
boundary conditions, or through the influence of (correlated) disorder, will therefore render 
the superconducting system effectively more three-dimensional. 

Uncorrelated point defects, e.g., oxygen vacancies in the copper-oxide planes, increase 
the fraction of entangled lines (hence the minus sign) by [Eq. ( |6.10|) 1 



n A = -An e 1 /4ir(k B Ty 



(6.98) 



clearly as a result of enhanced wandering of each flux line searching for an optimal path 
through the sample while trying to accomodate most profitably with the uncorrelated pin- 
ning centers. In very thick samples, the tilt modulus, however, remains unaffected by this 
kind of isotropically distributed disorder. The leading finite-thickness corrections to the 
fraction of pinned vortices remarkably display different power-law thickness dependences 
for periodic and open boundary conditions, respectively, see Eqs. (|6.14j) and fl6.15|) , 



An A (L) 



f An /A-nc\(k B TfL p.b.c. 
\ -9C(3)An /327reML 3 °- b - c - 



(6.99) 
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the corresponding expressions for the tilt modulus read 



(ni,)- 1 [1 - 3C(3)A/7re?c?L 3 ] p.b.c. 
(n^)- 1 [1 + 0(l/L 3 )] o.b.c. 



(6.100) 



For an ensemble of randomly tilted linear defects, i.e., nearly isotropic splay, we found that 
the disorder contribution to the fraction of entangled lines actually diverges [Eq. ( |6.38| )1, be- 
cause the thermal flux line wandering becomes superdiffusive as the vortices are "converted" 
along with the tilted defects. However, in thick samples the tilt modulus is unaffected by 
this isotropic disorder, and only in a finite sample does one get the following tilt modulus 
renormalization, 



n v — 1 
c 44 



{nhY 1 [1 - vrA/8e~ 2 cfL 2 ] 
(nei)" 1 



1 + vr 3 A/256e 2 cfL 2 



p.b.c. 
o.b.c. 



(6.101) 



The logarithmic divergence of the disorder-induced depletion TTX also occurs for parallel 
columnar defects (d± = 2), as long as the tilt angle a remains nonzero. (For this case, the 
divergence can presumably be removed by redefining the time-like direction to account for 
the vortex "drift".) In the spetial situation a = 0, Eq. ( 4.30 ) becomes 



nX = A£j/lQir 2 a(k B TY 



(6.102) 



and the corresponding finite-size corrections are nonzero only for open boundary conditions 
and read, according to ( |6.61| ), 



An A (L) = -(AG - ln2)An /27refc°L 



(6.103) 



with G ~ 0.916. For a tilt in the x direction with angle a > 0, the non-vanishing tilt 
modulus tensor components in a thick sample are [Eqs. ( |6.76| ), ( |6.77| )1 
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7r(/cB^) 2 (tana) 2 
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1 + e 2 (tana) /27rn a 

[1 + e 2 (tana) 2 /47rn a] 
1 



3/2 



7r(/c B T) 2 (tana) 2 \ [i + ^(t&n a) 2 /4im a 



1/2 



(6.104) 
(6.105) 



We remark again that this reduction of the tilt modulus by the correlated linear disorder can 
be reinterpreted as an increase of the mass anisotropy ratio mj_/m z ; the vortices are thus 
effectively rendered more three-dimensional by their interaction with the columnar defects, 
as should be expected. This mechanism increases the decoupling field in the presence of 
correlated disorder, above which the vortices change from lines, as they have been treated 
throughout this paper), to point-like "pancakes". For a — > 0, both the expressions (|6.104|) 
and ( |6.105| ) reduce to 



-■li 



nei 



Aij/87r 2 n a(k B Ty 



(6.106) 



which remains unchanged even for finite sample thickness in the case of periodic boundary 
conditions (toroidal geometry), while we find for open boundary conditions 
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-\L) 



mei 



7 ^-l 



1 - Ae 2 /87r 2 n a(fc B T) 2 + [16G + 21C(3)]A/87rl 2 c^ 2 



(6.107) 



For the infinitely thick sample, Eq. (|6.106| ) leads to an estimate for the Bose glass transition 
temperature from the high-temperature flux liquid phase, complementary to the estimate in 
Ref. We note that the average disorder potential for columnar defects is Vd = —Uob 2 /d 2 , 
where Uq ~ 0q/2(47tA) 2 denotes the defect potential well strength, b its width, and d the 
average distance between the columns [KJ. Accordingly, the defect correlator becomes 
Ao = U^b 4 / d 2 . One important modification has to be taken into account near the localization 
transition (more precisely, for temperatures above the depinning temperature T^ p ) , however, 
namely the effective thermal renormalization of both the disorder strength U (T) and the 
replacement of one factor of b 2 by the effective transverse localization length 1±{T), i.e.: 



A — > A(T) = U(Tyi ± (Tyb 2 /d 



(6.108) 



where 



U(T) = U (T*/T)' 



l ± (T) = d(T/T 



*\2 



(6.109) 



and the characteristic temperature is given by k^T* = b(iiUo) 1 / 2 ||40|| . Upon inserting these 
results, n Q = B/4> , and Eq. ( |6.93|) into Eq. ( |6.106| ), and collecting terms, we find that 



n v — 1 
c 44 



(T) 



1 - (Tb G /T) 4 



So C44 diverges as expected at the Bose glass transition temperature [^0 

Tbg = T*[0o/(47tA) 2 j B] 1 / 4 , 



(6.110) 



(6.111) 



in agreement with a similar analysis in Ref. f4H |. Although one should not trust the expo- 
nent for the vanishing of c^ -1 predicted by the approximate formula (|6.110|) , this analysis 
suffices to estimate the transition temperature. Upon generalizing to tilted linear disor- 
der, and utilizing Eqs. ( |6.104| ) and ( |6.105| ) we can locate the instabilities for tilt in the x 
and y directions, respectively, as a function of the tilt angle a, or instead introducing the 
dimensionless parameter 



u(a) = e 1 (tana) /Aimoa 



and find 



T x (a) 


4 _ 2 


Tbg 


u(a) 


T y (a)' 


4 2 


Tbg 


u(a) 



2u(a) 



- 1 



[l + «(a)]V2 



(6.112) 

(6.113) 
(6.114) 



The temperatures T x and T y , where the instabilities occur, are plotted in Fig. |]as a function 
of the parameter u. Notice that while always T y (a) > 0, T x (a c ) = for u c — (1 + v / 5)/2, 
suggesting that the Bose glass remains stable towards tilt up to this critical tilt angle a c . 
For a < a c , the system will thus not respond to transverse magnetic fields, i.e., display 
a transverse Meissner effect |4~0" |. The curve T x (a) can therefore be viewed as describing 
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the phase boundary (localization transition line) in a phase diagram with the thermody- 
namic variables temperature and transverse magnetic field. However, within the Bogoliubov 
(Gaussian) approximation used here, we do not find the cusp in T x (a) that is predicted by 
the Bose glass scaling theory [[|(| . 

For the case of two families of symmetrically tilted (angles ±a) columnar defects, the 
fraction of pinned lines remains finite because the leading logarithmic divergences for both 
directions precisely cancel; Eq. ( |6.85| ) yields 



1/2 



nK = Ael/16n 2 a(k B T) 2 1 + e 2 (tana) 2 / '4irn a 

and the tilt modulus components (for L — > oo) deriving from Eqs. ( |6.89|) and (|6.90| ) are 

1 



(6.115) 
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(6.116) 
(6.117) 



As before, we can investigate the instabilities of this system towards external tilt, and find 



[T x (a)/T BG ] 4 =[l+u(a)}- 3 / 2 
[T»/T BG ] 4 = [! + «(«)] 



-i/2 . 



(6.118) 
(6.119) 



note that both T x (a) are always T y (a) positive (see Fig. ||), and the Bose glass remains 
stable irrespective of the tilt angle a. 

Turning to parallel defect planes (c?j_ = 1), oriented in the yz plane (i.e., a = 0), the 
reduction in the boson order parameter ( |4.38| ) is 



WX = Aii/QA^ 2 nl /2 a z ' 2 k B T . 



(6.120) 



For any nonzero tilt angle, we again find a divergence, as in the case of linear defects. 
According to Eq. ( |6.60|) , the leading finite-L corrections to (|6.120|) for a = are in the case 
of open boundary conditions 



An A (L) = -(tt - l)An /2'Ki 2 l c\L 



(6.121) 



while no such finite-thickness contributions occur if periodic boundary conditions are ap- 
plied. Tilting the defect planes in the x direction with angle a, the components of the tilt 
modulus [see Eq. |6.75|| read 
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16(7rn af/ 2 k B T [l + gf(t ana) 2 /47moa 



(6.122) 
(6.123) 



i.e., the tilt modulus in the y direction transverse to both the tilt and the magnetic field 
remains unaffected by the disorder. When a — > 0, Eq. ( |6.122| ) becomes 



->u —1 



(ni)- 1 (l - Ail/lQ(nn a) 3 / 2 k B T 



(6.124) 
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only for open boundary conditions will there be corrections for finite sample thickness, which 
we have worked out in Eq. fl6.66j ) to be 

c^-^L) = (ni)- 1 [l - Ail/16(7rn a)V 2 k B T + (vr + A)/Al\c\L\ . (6.125) 

As for the linear defects, we can estimate a critical temperature from Eq. ( |6.124j ). The 
"bare" disorder correlator for planar defects is U^b 2 /d, where b denotes the diameter of 
the defect potential, and d the average defect distance. The thermal renormalizations for 
T > T dp now read |40] 

A — ► A(T) = U(T) 2 b 2 /d , (6.126) 
U{T) = U (b/d) 2/3 (T*/T) 2/3 , l ± (T) = d(b/d) 2/3 (T/T*) 4/3 , (6.127) 

implying that A(T) = U^b 3 /d 2 , independent of T. This yields 



cl x ~\T) = {nhV [1 " Tbg/T\ , (6.128) 

and the instabilitiy occurs at the Bose glass transition temperature marking the onset of 
localization to the defect planes 

f BG = T*(b/Ad) 2 ((P /2n\ 2 B) 3/2 . (6.129) 

As a function of tilt angle a, or equivalently, the parameter u introduced in Eq. ( |6.112j ), we 
furthermore find 

T ° - 1 - 2 " (a) (6.130) 



f BG [l + u{c)fr- ' 

see Fig. [| In this approximation, the critical angle beyond which the Bose glass phase of 
localized vortices becomes unstable towards tilt is given by u c = 1/2. 

Upon generalizing to an ensemble of two families of symmetrically tilted planar defects 
(in thick samples), the divergences present for each of the families with positive and negative 
tilt angle, respectively, again cancel, and the resulting average density of pinned flux lines 

(UD ^ 

n£ = Ae^/647r 3/2 no /2 'a 3/2 k B T [l + e^(tan a) 2 / '47m aY ^ , (6.131) 

and the tilt modulus [from Eq. (|6.88| )1 along the tilt direction reads 

^T~T_ 1 (i Ae ~i l-6~f(tana) 2 /87m a \ 

xx " nh \ 16(nn a)^k B T [1 + gf(tana)V47rn a] 5/ V ' 1 ' ' 



while of course c v yy 1 = (^ei) 1 . The angle dependence of the instability in c v xx 1 is readily 
found to be 

T x _ l-u(a)/2 

f BG -[l + «(a)F' [bAdd) 

compare Eq. ( |6.130| ). As is shown in Fig. [], in a system with two families of symmetrically 
tilted planar defects the stability region of the Bose glass in the phase diagram is enhanced, 
the critical angle now being determined by u c = 2. 
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VII. SUMMARY AND DISCUSSION 



In this paper, we have investigated the influence of weak uncorrelated and correlated 
disorder on superfluid bosons as well as magnetic flux liquids in high-T c superconductors, 
exploiting the formal analogy of the statistical mechanics of directed lines in d+1 dimensions 
with the path integral representation of the quantum mechanics of d-dimensional particles. 
We have specifically addressed the important issue of how different boundary conditions in 
the directed-polymer case modify the results in systems with finite thickness L. 

For disordered boson superfluids, we have generalized the previous studies by Huang 
and Meng JT7|] , and Giorgini, Pitaevskii, and Stringari on the effect of point disorder 
on the depletion as well as the superfluid density to the case of correlated disorder, i.e., 
linear or planar defects. In Sec. |TV| , we have derived general results for the density, current, 
and vorticity correlation functions, as well as the condensate fraction and the normal-fluid 
density. Explicit formulas for the contributions by pointlike, linear, and planar disorder, 
as well as randomly positioned and oriented line defects, to the depletion and normal-fluid 
density in two and three dimensions are summarized in Sec. [IV Q . We have also found that 
these results for static disorder are independent of temperature (at least to first order in the 
defect correlator). 

In the case of magnetic flux line liquids in high-temperature superconductors, we saw 
that different boundary conditions can actually change the effects of disorder qualitatively. 
Even in the absence of disorder, we found that a homogeneous condensate fraction per- 
sists at finite sample thicknesses if open boundary conditions are applied, contrary to the 
case of periodic ("bosonic") boundary conditions, where long-range order is destroyed by 
"thermal" fluctuations in this (2 + l)-dimensional system. We have physically related the 
qualitatively and quantitatively different results for periodic and open boundary conditions 
to the enhanced thermal wandering of flux lines near the sample boundaries in the latter 
case, see Fig. ||. In Sec. [V|, we have studied the situation for open boundary conditions in 
some detail, looking at precisely those quantities which are the analogs to those studied for 
the Bose superfluid, namely density and tilt correlation functions, the "boson order param- 
eter" , which is measure of the fraction of entangled lines, and the tilt modulus. In Sec. [VT|, 
we have compared our findings in two and three dimensions for the two kinds of bound- 
ary conditions we explored, specializing to (in the flux line language) point defects, nearly 
isotropically splayed linear disorder, parallel untilted and tilted extended defects (lines and 
planes), as well as two families of symmetrically tilted extended defects. We have thus sig- 
nificantly generalized earlier work on weak point disorder in infinitely thick systems [|51|] in 
two important respects, namely by (i) studying correlated disorder as well, and (ii) taking 
carefully into account the role of the more realistic open boundary conditions. 

In Sec. [VI F| , we provide explicit results for the average fraction of entangled flux lines and 
the tilt modulus, and elucidate some of their experimental consequences. E.g., we have com- 
mented on how the defect influence on the tilt modulus effected by anisotropic (correlated) 
disorder may be reinterpreted as a renormalization of the effective mass anisotropy, illustrat- 
ing how columnar and planar defects render the system effectively more three-dimensional; 
this renormalization clearly increases the decoupling field, above which the vortices degen- 
erate from three-dimensional objects (lines) to weakly coupled two-dimensional ones (pan- 
cakes). Furthermore we have discussed how the renormalized tilt modulus may be used to 
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estimate the location of the instability towards a localized phase from the high-temperature 
side of the phase diagram, and to study the stability of this Bose glass with respect to 
transverse magnetic fields ("transverse Meissner effect"). 

Our entire investigation was based on the Gaussian approximation, and nonlinear effects 
have been neglected. Analysis of these nonlinearities, as well as disorder which couples 
directly to the tilt field (corresponding to a boson current) are interesting avenues for future 
investigation. We have restricted our study to local (in z) and isotropic repulsive vortex 
interactions, which is a fair approximation for low magnetic fields. It would be interesting 
to extend our results to the high-field regime, perhaps along the lines suggested in the 
Introduction. 
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APPENDIX A: DERIVATION OF THE DISORDER EFFECTS ON THE TILT 
MODULUS USING AFFINE TRANSFORMATIONS 



In this first Appendix, we provide yet another derivation of the disorder-induced renor- 



malization of the tilt modulus of a flux line liquid, utilizing an affine transformation [44 



which is applicable in the (thermodynamic) limit of thick samples L — > oo, and for disordered 
systems in which translational invariance is statistically restored. Upon applying a uniform 
tilt ( |2.62| ), the transformed action in polar representation reads [see Eqs. ( p.55| ),( |236|) 1 



Sy, 0'] = S[ir, 0] + J dzj d d r 



■ fc B T 
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h-V6(r,2) 



~2 



6Vn(r,z) 



(n + vr(r, z) 



(Al) 



where S[n, 0] denotes the pure part (in the flux line notation) of the action (|3.2| ). This can 
also be written as 



sy,e'\ = sy,e'\ - idz' id d r' 



5V D (r'-hz'/e 1 ,z') + ^ 



(n + 7r'(r', z')) . (A2) 



The partition function for the tilted system therefore becomes (we henceforth omit the 
prime) 



In Z gr = In Z gr + In ^ exp 



dz / d d r 



5V D (r - hz/h,z) 
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(A3) 



where the thermodynamic average (...) is to be performed using the pure action, and the 
vortex contribution to the tilt modulus [see Eq. can then be inferred from 



k B T 5 2 lnZ. 



n 2 LVt 5hi5hj 



(A4) 



h=0 



In order to compute the thermodynamic and quenched disorder averages, we perform a 
double cumulant expansion, which yields 



In Z RI = In Z + - , - h 2 + 



hk B T 2(k B T) 2 J (27r) d J 2tt 



d d q f dq 



A(q,&;-q, ~q z )h+ (^z~ 



(A5) 



here, S(q,q z ) denotes the density correlation function of the pure system ( [3.151 ). Assuming 
now that the quenched disorder average restores translational invariance statistically, 

A(q, q z ] q\ q' z ) = A(q, q z ) fi5 q ,_ q . LSq,,-^ , (A6) 
and expanding the disorder correlator to second order in the tilt, one obtains with Eq. (|A4[) 
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If the system is statistically isotropic (with respect to the d transverse dimensions) as well, 
Eq. ( |A~8f) reduces to J|4| 
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Inserting 
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we finally arrive at 
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compare Eq. ( |4.15|) in Sec. [IV A . 

Eq. (|A7|) immediately shows that uncorrelated point defects, A(q, q z ) = A, and "almost 
isotropically" splayed columnar defects, A(q, q z ) = A/q [[H], have no influence on the tilt 
modulus in thick samples at all, which is true for any disorder whose correlator is independent 
of q z . For parallel extended defects, tilted in the x direction, 

A(q||,q_L,<? 2 ) = AVL\\5^ L5 q ^ qx 

tan a > 

(A12) 

one sees that the disorder does not change the tilt modulus along the defect directions, 

c| = ne x , (A13) 
while in the perpendicular directions, compare Eq. fl6.74| ), 
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As discussed in Sec. [VI F| , a divergence of the transverse tilt modulus in a certain direction, 
0, marks a localization transition to a (generalized) Bose glass phase; Eq. ( |A14| ) 



„v — 1 



can thus be employed to obtain a "high-temperature" estimate for the Bose glass transition 
temperature [|4*4" |. 
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APPENDIX B: PHASE ONLY APPROXIMATION 



Considering the two types of fluctuations in our system, described by the harmonic action 
(P-2|), it may be tempting to integrate out the density fluctuations, which acquire a mass 
Vo, see Eqs. (|4.1|) and (|3.8| ), and only retain the massless phase modes; certainly this should 
correctly describe the leading terms in a long-wavelength expansion. In this Appendix, 
we study this simple phase-only approximation (with periodic boundary conditions), and 
explain its deficiencies. 

We thus define the effective action S^frfO] by 

e - Seff [0|A = |p Me -So[,,ep (B1) 
The purely quadratic functional integral is readily evaluated, with the result 

S cS [&] = COnst. + £ p(g)2 [^m + e B(g)V^ 2 ] (-q ) -^m)0(q,^m) 

~T~ — 2hu m Q(-q, -uj m ) + SV D (-q, -u m ) 5V D (q,u m ) \ ; (B2) 
2m e B {q) 2 L J I 



shifting the B field as in Eq. (|4.3|) , one then finds for the phase correlation function precisely 
the result (|4.4j) . However, as density fluctuations have now been neglected entirely, the 
Green's function is simply 



G(q, u m ; q', urn') « n (O(q, u m )Q(q', u m >)) , (B3) 
and similarly the current correlation function is merely given by the first term in Eq. ( |3.21 ), 
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where, in order to be consistent with our long-wavelength expansion, we should use €b(?) ~ 
hciq. Thus all the nonlinear terms, i.e., the vorticity contributions fl3.27| ) present in the 
"full" theory ( |4.1|) have disappeared, and as the current correlation function is therefore 
purely longitudinal, see Eq. ( |3.28|) , the normal-fluid density will necessarily be zero at all 



temperatures in the phase-only approximation, 



Pn{T) = . (B5) 

Using Eq. ( |3.51| ), we can also compute the average depletion, 

_ mci f d d q If 2 

n-n = — / — — - 1 1 



2 J (27c) d q V - 1 

f d d q 1 LOmUm' A(q, u m ; -q, u m >) , Rfi x 

+ h 2 n J {2«Y {(3hf ±, (oft + clq 2 )^ + c\q 2 ) ' 1 J 
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For the pure system, this yields the at least qualitatively correct result that the zero- 
temperature depletion diverges logarithmically in one dimension; however, for d > 2, 
n(T = 0) cannot be determined from Eq. (|B6j) . On the other hand, the finite-temperature 
corrections An(T) turn out to be identical to the leading contributions ( |3.58| ) obtained for 
the "full" theory in the phonon approximation [i.e., neglect any terms of order q 2 in the 
integrand of Eq. ( [OS )]. In the same manner, the disorder contribution to ( p6|) yields the 



leading terms in a long- wavelength expansion for the depletion, whenever these do not van- 
ish. E.g., for uncorrelated ("point") disorder in space and imaginary time, one finds [see 

n f d d q coth[/3e B (g)/2] / (3e B (q) 



and for the case of "almost isotropic splay" accordingly, cf. Eq. (|6.35| ), 

— — - n / d d q coth[/3e B (g)/2] ( (3e B {q) \ 

= -^iW? ieM V ~ sinh[/? eB (g)]J ' (B8) 

while for "tilted extended defects" the result is temperature-independent [see Eq. (|6.72|) 1, 



_ A ^o f d d q ± 

UA " h 2 J (2tt)^ [e B (q ± ) 2 /h 2 + q 2 v 2 ] 2 ' 1 ' 

All of these expressions correctly represent the leading contributions of the results obtained 
in the "full" theory where the density fluctuations are taken into as well. However, at 
T = the integrals ( |B7| ) and ( |B8| ) become ill-defined when the phonon approximation is 
employed, as is required for the sake of consistency here. Note in addition that for the case of 
"untilted" correlated disorder (along the r direction), and for two families of "symmetrically 
tilted extended defects" the lowest- wavenumber contributions actually cancel, leaving the 
very poor result n&(T) = 0. 

It is actually this very mechanism which leads to a vanishing normal-fluid density in 
the phase-only approximation, namely a cancellation of the leading-order terms in the 
long-wavelength expansion. This can be made explicit by actually calculating the aver- 
age momentum density in a system with walls moving relative to the superfluid; a quick 
computation shows that even with disorder 



(g(r,<r)> = O + 0[(q-v)1 , (BIO) 

quite in accord with the fact that within the framework of the phase-only approximation 
there is no vorticity in the system, and hence p n {T) = for all temperatures. 

Summarizing our experiences, we have realized that the phase-only approximation ( |B2|) 
is sufficient only for the evaluation of the longitudinal current correlations, and typically also 
for the leading finite-temperature corrections to the condensate fraction. It produces quite 
unsatisfactory results for the transverse current and vorticity correlations, however, which 
are crucially affected by their coupling to density fluctuations, and hence fails to detect any 
nonzero normal-fluid density at finite temperatures. 
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APPENDIX C: OPEN BOUNDARY CONDITIONS: CORRELATION 
FUNCTIONS AND DISORDER CONTRIBUTIONS 



In this Appendix, we derive the density and phase correlation functions in Gaussian 
approximation for a disordered flux line system described by an action of the form ( |3.2|) , 
with open boundary conditions. This calculation follows and generalizes Refs. [[31]] and ||51|| . 
i.e., the disorder contributions will be treated to first order in the disorder correlator A, 
and the open-boundary correlators will be expressed in terms of matrix products of the 
corresponding correlation functions obtained with periodic boundary conditions. Although 
the physical application we have in mind is that for magnetic flux liquids, we retain the 
boson notation here in order to facilitate comparison between the two cases. 

Our starting point is Eq. ( |3.72|) , where the boundary conditions (|3.70|) have been explic- 
itly taken into account through the insertion of a static constraint field A(r). Introducing 
the vector 

«(q,a^)=( (%( ° Wm) )+<nA(q) ) (CI) 

we can write the new effective harmonic action as 
^o[A, T] = ^ [T T (-q, -u; m )A(q,^ m )T(q,w m ) -2T T (-q, -cj m )M(q,cj m )] , (C2) 



where A(q, u m ) is the matrix Q3.8|), obeying A (— q, — u m ) = A(q, tu m ). Transforming (with 
Jacobian 1) to new fields T(q, u m ) = T(q, ui m ) — A _1 (q, u m )u(q, ui m ) [see Eq. ( |4.3|) 1 , the 
effective action So [A, T] is purely quadratic in T, and hence the functional integral over these 
variables is readily performed, leading to 

(T(q, w m )T T (q / , u m t))= hA^(q, u n ) f3hVt 5 q ,_ q /<5 m - m ' 

+A _1 (q, uj m ) (u(q, u m )u T (q', LU m >))x A _1 (-q', -u m t) , (C3) 

with the inverse matrix A _1 (q, u m ) given in Eq. ( |3.9| ), and 

(t*(q, u m )u (q , Wmf )> A = JV[X(q)}e-^im ' (C4) 

where 

Sy[\] = u T (-q,-u m )A~\q,u m )u(q,u m ) . (C5) 

At this point we caution the reader to note that Eq. (|C4|) does not simply reduce to the 
negative of the first line in Eq. ( |C3| ) , for the functional integral is over the static field A(q) 
only, and not over the dynamic variable u(q, ou m ); the fluctuations will only vanish at the 
special points r = and r = (3%. 

We thus have to evaluate Eq. (P4|) very carefully. Upon defining 



A-\q) = A^q, r = 0) = -L £ A-^q, u m ) , (C6) 

pn m 
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the inverse of which is explicitly written down in Eq. ( P-74Q , and furthermore 



(q) = ^? A " 1(q '" m) ( 





SV D (q, uj m ) 



(C7) 



and then again shifting (with Jacobian 1) the integration variable A(q) = A(q) — 
iA(q)a(q)/h, it becomes a straightforward task to perform the purely quadratic path inte- 
gral over A, with the result 

(u(q, u m )u T (q', u m >)}\ = -ftA(q)Q6 q _ q # + 5V D (q, uJ m )5V D (q', u m ,) ( ® J 
-A(q)-^-^A- 1 (q,^ n )5V D (q,w n ) ^ J ^ 5Vb(q', Wm') 
-5V D (q,u; n ) (° J J -1 £ A^-q', -u>„0*Vb(q', ^)A(q') 

+ A (q)^-E A_1 (q^n)^(q^n) (o ?) ^EA^-^-^w^OAtq') . (cs) 

Finally performing the quenched-disorder average, and then collecting the various terms, 
we find for the matrix of the averaged density and phase correlation functions 



(T(q,^ m )T T (q',^ m /)) = (2ir) d 5(q + q')hA \q,uj m ) (3M m _ m i - A(q)A \q,-w. 



+A 1 (q,u; ri 



Z)|[^m,n- A(q)A X (q,< 







o 





X 



X 



Phfim'n' ~ A(q')A ^q', LU n >)] A(q, u; n ; q', u n >) }A 1 (-q / , -w n 



(27r) d 5(q + q') [/?M m ,_ m / - A x (q, w TO )A(q)] TiA x (q, -w„ 
— - [Ph$m,n ~ A" 1 (q,^ m )A(q)] A -1 (q,av) ( J ) X 



xj [3h5 m/ y - A _1 (q',^ m /)A(q') A -1 (q',W n /) J J A(q, w n ; q', w n > 
It can be easily checked that indeed the fluctuations vanish at the boundaries, 



(C9) 



(T(q, r = 0)T^(q', r = 0))= (T(q, r = ^)T^(q', r = /%)} 

' J2 (T(q,c m )T^(q',^))=0 



(CIO) 



as required by Eqs. ( |3.70| ). Also, if all the second terms in the square brackets, which break 
translational invariance in r, are disregarded, we recover the two-point correlation functions 
34|)-( |4.6|) for periodic boundary conditions obtained in Sec. [IV A| . 
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APPENDIX D: MATSUBARA FREQUENCY SUMS AND MOMENTUM 

INTEGRALS 



For the reader's convenience, we finally provide a list of useful explicit formulas for the 
(bosonic) Matsubara frequency sums required in the above calculations, as well as some 
definite integrals needed for the momentum integrations. 



1. Bosonic Matsubara frequency sums 

The following fundamental formula for bosonic and fermionic Matsubara frequency sums 
(n = 0, ±1, ±2, . . .) is readily obtained by performing an appropriate contour integration in 
the complex plane (for details see Ref. [0], p. 248 ff.) 

1 v e lU)nT _ _ e Te/h i bosons with uj n = 2rm / (3% , 
T - Ph^iun-e/h ~ T e^ T l tOT j fermions with u n = {2n + 1)tt / /3% . {D1} 

Consequently, 

J_ v e iuJ " T h_ e V- T M e ± e Te / h h_ j coth (/3e/2) (bosons) 

j3h V a* + e 2 /Ti 2 ~ 2e T 1 ~~ ' 2e j tanh (/3e/2) (fermions) T ~* ' ( j 

Now specializing to bosons, i.e., periodic boundary conditions in imaginary time and 
therefore Matsubara frequencies with even integers, we may write down the following se- 
quence of formulas (for r > 0): 



(3h „ iu n + e/k (3Ti ^ iu n — e/% 



(3h 5 iu n - e/h f5h ^ iu n + e/% - 1 ' 



1 v e-^" r ftcosh[(/%-2r)e/2ft] 

V ^ 2 + e 2 /fr 2 ~ 2esinh (/3e/2) ' ' °' 

J_ ^ zcu w e-^" T _ sinh [(/ffi - 2r)e/2ft] , „ 

/3ft^^ + e 2 A 2 2sinh(/3e/2) 



One then immediately sees that (m, n = 0, ±1, ±2, . . .) 

1 ^ m -^)e-^-+^ 

and hence 

1 (uj m m n + e 2/h 2 )e-^ + ^ eg _ 1 

(/3ft) 2 ^ K + e 2 /ft 2 )(^ 2 + e 2 /ft 2 ) (e* - l) 2 4 [sinh (/?e/2)] 2 ' 1 J 

independent of r. 

Differentiation of Eqs. (|D5|) and ( p6|) with respect to e furthermore yields 
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E 



n 3 



ph V K + e 2 /^ 2 ) 2 4e 3 sinh (/3e/2) 



cosh 



{ph - 2r)e 



— — sinh 

n 



2h 

(ph - 2r)e 



E 



0/i V K + e 2 /^ 2 ) 2 4e sinh (pe/2) 
and combining these with Eq. (03) gives 



r cosh 



2h 

(]3h - 2r)e 
2h 



(3e cosh (re/h) 
+ Y sinh (/3e/2) 
/3ft sinh (re/ft) | 
T sinh (/3e/2) J ' 



1 ^ uj 2 n e~ iu > nT h 



Ph^ {ul + ey?i 2 ) 2 4esinh(/3e/2) 



cosh 



1 (uj 2 n -e 2 /h 2 )e- 1 ^ 
Ph\ (^ + e 2 /ft 2 ) 2 



2 sinh (/3e/2) 



re • i 
+ — sinh 

n 



t sinh 



(Ph - 2r)e 
2fr 

(Ph-2r)e 

2h 

(ph - 2r)e 
2~h 



Pe cosh (re/h) 
Y sinh (Pe/2) 
Ph cosh (re/h) 
y sinh (/3e/2) 



(D9) 
(D10) 



(Dll) 
(D12) 



2. Momentum integrals 

We finally list some useful formulas for some of the momentum integrals required for 
Bose systems in d dimensions, see Refs. [0,^]: 



In 2 for d = 2 , 

l-2-( d - 2 Ar(d-l)((d-l) for d>2, 



-dx= T(d-l)C(d- 1) , 
o &- — 1 

x d-2 

dx 

o e 11 + 1 

^-2- cte= 2~ (d ~ 2) r(d) C(d - 1) , 
o sinh x 

rfx= 2-( d ~ 2 ) ( 1 - 2-^ 2 )) Y(d) ((d - 1) 
o cosh x v ' 

x d ~' 



o sinh x 

rj.d—1 



dx= 2(1- 2- d ) Y(d) ((d) 
dx=2Y(d)i(d) . 



(D13) 
(D14) 

(D15) 
(D16) 
(D17) 
(D18) 



lo coshx 

Here, ((n) = E£Li 1/*" and £(n) = ET=o(- 1 ) h I \ 2k + l T- We note that the C function at 
even arguments is related to Bernoulli's numbers, 



C(2*) 



2fco2A;-l 



7l 2k 2 



(2k)\ 



B k , C(0) = 4,C(2) = ^-,C(4) = ^,C(6) = ^ 



while the £ function at odd arguments is connected with Euler's numbers, 



C{2k + 1) 



IT 



2k+l 



E k , e(i) = f,e(3) = ^... 



4 fe + 1 (2fc) 

and C(3) « 1.202, £(2) = G « 0.915965 (Catalan's constant || 



(D19) 



(D20) 
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TABLES 



TABLE I. Superfluid boson - flux liquid analogy. 



Superfluid bosons m % f3h V{r) n [i iv p s 

Vortex liquid l x T L 2^K^FJ\) B/fo (H - H^q/Att h/i[ p 2 c~ AA 



TABLE II. Varieties of disorder. 



Vortex liquids in three dimensions 
Point disorder (e.g., oxygen vacancies) 
Columnar disorder (e.g., screw dislocations, 

columnar damage tracks) 
Planar disorder (e.g., twin boundaries, 

grain boundaries) 



Superfluid Helium films 

Space- and time-dependent random potential 
Amorphous substrate modelled by uncorrelated 

"point-like" random potential 
"Smooth" or crystalline substrate with line 

disorder deposited microlithographically 
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FIGURES 



FIG. |l|. Schematic representation of the Abrikosov flux lattice (a), a disentangled flux 
liquid (b), and an entangled flux liquid (c). In the boson picture, the directed lines map 
onto the imaginary-time particle world lines, and these phases correspond to a crystalline 
solid (a), a normal liquid (b), and a superfluid (c), respectively. 

FIG. |2|. Schematic plot of the time-averaged cross section of a thermally fluctuating flux 
line in a slab with free boundary conditions. The bulk part away from the boundaries has a 
width determined by the square of the bosonic "wave function" , which is a Gaussian for a free 
random walker. Near the top and bottom surfaces, thermal wandering is enhanced, reflected 
by a probability distribution determined just by the "wave function" (not its square). 

FIG. |3|. Lowest-energy solid phase contribution to the correlation function G(r, z; r', z'), 
which inserts a flux head and tail into a crystalline vortex array. Dashed lines represent a 
row of vortices slightly behind the plane of the page. In (a), a vacancy is created at "time" 
z, which then propagates and is destroyed at "time" z'. Interstitial propagation from z' to z 
is shown in (b). The energy of the "string" defect connecting the head to the tail increases 
linearly with the separation in both cases and leads to the exponential decay of G(r, z; r', z'). 

FIG. £|. Critical temperatures T x (a)/T BG (solid line) and T y (a)/T BG (dashed) where the 
tilt moduli c v xx ~ x and c yy ~ l in a system with parallel columnar defects tilted along the x direc- 
tion vanish, respectively, as a function of the dimensionless parameter u = e\(tan a) 2 / Aimoa. 
The curve T x (a) /Tbg may be viewed as an estimate for the phase boundary of the Bose 
glass, which remains stable towards external tilt as long as a < a c (transverse Meissner 
effect). 

FIG. ||. Tilt instability temperatures T x (a)/T B Q (thick solid line) and T y (a)/T BG 
(dashed) for a system with two families of symmetrically tilted columnar defects. For 
comparison, T x (a) /Tbg (thin solid line) for a single family of tilted linear disorder is also 
depicted. 

FIG. |6|. Tilt instability temperature T x (a)/TBG for both parallel planar defects uni- 
formly tilted along the x direction (solid line) and two families of symmetrically tilted planes 
(dashed). The critical tilt angle marking the instability of the Bose glass phase of localized 
flux lines is shifted upwards in the symmetric-splay case. 
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- FIG. 7 in Ref. 
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FIG. 1. Schematic representation of the Abrikosov flux lattice (a), a disentangled flux liquid 
(b), and an entangled flux liquid (c). In the boson picture, the directed lines map onto the 
imaginary-time particle world lines, and these phases correspond to a crystalline solid (a), a normal 
liquid (b), and a superfluid (c), respectively. 




FIG. 2. Schematic plot of the time-averaged cross section of a thermally fluctuating flux line 
in a slab with free boundary conditions. The bulk part away from the boundaries has a width 
determined by the square of the bosonic "wave function", which is a Gaussian for a free random 
walker. Near the top and bottom surfaces, thermal wandering is enhanced, reflected by a probability 
distribution determined just by the "wave function" (not its square). 
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- FIG. 12 in Ref. 
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FIG. 3. Lowest-energy solid phase contribution to the correlation function G(r, z; r', z'), which 
inserts a flux head and tail into a crystalline vortex array. Dashed lines represent a row of vortices 
slightly behind the plane of the page. In (a), a vacancy is created at "time" z, which then propagates 
and is destroyed at "time" z' . Interstitial propagation from z' to z is shown in (b). The energy of 
the "string" defect connecting the head to the tail increases linearly with the separation in both 
cases and leads to the exponential decay of G(r, z\ r', z'). 




FIG. 4. Critical temperatures T x {a)/T-QQ, (solid line) and T y {a) /T-qq (dashed) where the tilt 
moduli c v xx ~ x and c v yy ~ x in a system with parallel columnar defects tilted along the x direction 
vanish, respectively, as a function of the dimensionless parameter u = ef (tan a) 2 /Annoa. The 
curve T x (a) /Tbg may be viewed as an estimate for the phase boundary of the Bose glass, which 
remains stable towards external tilt as long as a < a c (transverse Meissner effect). 
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T,(a) 



BG 1.0 




FIG. 5. Tilt instability temperatures T x {q)/T^q (thick solid line) and T y {a)/T-QQ (dashed) for 
a system with two families of symmetrically tilted columnar defects. For comparison, T x (a) /Tbg 
(thin solid line) for a single family of tilted linear disorder is also depicted. 
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1.2 



T»(a) 




FIG. 6. Tilt instability temperature T x {a) /Tbg for both parallel planar defects uniformly tilted 
along the x direction (solid line) and two families of symmetrically tilted planes (dashed). The 
critical tilt angle marking the instability of the Bose glass phase of localized flux lines is shifted 
upwards in the symmetric-splay case. 
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